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ABSTRACT 

We present the power spectrum of the reconstructed halo density field derived from a sample 
of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Re- 
lease (DR7). The halo power spectrum has a direct connection to the underlying dark matter 
power for k ^ 0.2 h Mpc -1 , well into the quasi-linear regime. This enables us to use a fac- 
tor of ~ 8 more modes in the cosmological analysis than an analysis with k max = 0.1 h 
Mpc" 1 , as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 
2006). The observed halo power spectrum for 0.02 < k < 0.2 h Mpc 1 is well-fit by our 
model: \ 2 = 39.6 for 40 degrees of freedom for the best-fitting ACDM model. We find 
f2 m /i 2 (n s /0.96) 013 = 0.141~^q q°2 f° r a power law primordial power spectrum with spectral 
index n s and Vl^h 2 = 0.02265 fixed, consistent with CMB measurements. The halo power 
spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to 
an effective distance to z = 0.35: r s /D v (0.35) = 0.1097t°'.oo42- Combining the halo power 
spectrum measurement with the WMAP 5 year results, for the flat ACDM model we find 
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f2 m = 0.289 ±0.019 and H = 69.4±1.6kms 1 Mpc 1 . Allowing for massive neutrinos in 
ACDM, we find m v < 0.62 eV at the 95% confidence level. If we instead consider the ef- 
fective number of re lativistic species N e f f as a free parameter, we find N e f / = 4.8^7. Com- 
bining also with the Kowalski et alj d2008f) supernova sample, we find f2 tot = 1.011 ± 0.009 
and w — —0.99 ±0.11 for an open cosmology with constant dark energy equation of state 
w. The power spectrum and a module to calculate the likelihoods is publicly available at 
|http : //lambda . gsf c . nasa . gov/toolbox/ lrgdr/| 

Key words: cosmology: observations, large-scale structure of Universe, galaxies: haloes, 
statistics 



1 INTRODUCTION 

The past decade has seen a dramatic increase in the quantity and 
quality of cosmological data, f rom the discovery of cosmological 
accel eration using supernovae dRiess et aT]| 19981 ; IPerlmutter et al.l 
1999) to the precise mapping of the cosmic microwave back- 
round (CMB) w i th the Wilkinson Microwave Anisotropy Probe 
Page et all 120031 ; iNolta et all 120091) to the detection of the im- 
print of baryon acoustic oscillations (BAO) i n the early uni - 
verse on galaxy clustering ( lEisensteinetalj2 005; Cole et al.l2005h . 
Combining the most recent of these three cosmological probes, 
iKomatsu et al.l d2009h detect no significant deviation from the min- 
imal flat ACDM cosmological model with adiabatic, power law 
primordial fluctuations, and constrain that model's parameters to 
within a few percent. 

The broad shape of the power spectrum of density fluctuations 
in the evolved universe provides a probe of cosmological parame- 
ters that is highly complementary to the CMB and to probes of 
the expansion history (e.g., supernovae, BAO). The last decade has 
also seen a dram atic increase in the sc ope of galaxy redshift sur- 
veys. The PSCz dSaunders et ai1l2000h contains ~ 15000 IRAS 
galaxies out to z = 0.1, th e 2dF Galaxy Redshift Survey (2dFGRS; 
IColless et al]|200ll . l2003l) collected 221,414 galaxy redshifts with 
median redshift 0.11, and the Sloan Digital Sky Survey (SDSS; 
lYorketalJl2000l) is n ow complete with 929,555 galaxy spectra 
( Abazajian et al. 2009) including both main galaxies ({z) ~ 0.1; 
Strauss et al. 2002) and Luminous Red Galaxies (LRGs; z ~ 0.35; 
Eisenstein et alj |2001). To harness the improvement in statistical 



power available now from these surveys requires stringent under- 
standing of modeling uncertainties. The three major components 
of this uncertainty are the n on-linear gravitational evolution of 
the matter density field (e.g.. IZel'dovichlll970l ; iDavis et alj|l977l ; 
iDavis & Peebles! Il977l) . the relationship between th e galaxy and 
underlying matter density fields ( "galaxy bias", e.g.. Kaisej l 19841 ; 



Rees 1985; Cole & Kaiser 1989), and redshift space distortions 
(e.g.. iKaiseill 19871 : IDavis & Peebles! 1 19831 and iHamiltonll 19981 for 
a review). 

Several major advances have enabled previous analyses of 
2dFGRS and SDSS to begin to ad dress these complic ations. 
Progress in iV-body si mulations (e.g.. |Heitmann et alj|2008l) . an- 
alytical methods (see ICarlson et al j 120091 for an overview and 
com parison of many recent methods), and co mbinations thereof 
(e.g. JSmith et al, 2003; Eis enstein et alj|2007b[) have allowed sig- 
nificant progress in the study of the non-linear real space mat- 
ter power spectrum. Recent power spectrum analyses have ac- 
counted for the luminosity dependence of a sca le independent 
galaxy bias dTegmark et al.l2004al ; lcole et alj|2005l) . which can in- 
troduce an artificial tilt in P(k) in surveys which are not volume- 
limited dPercival et alj2004l) . ICresswell & Percivall d2009l) have re- 



cently examined the scal e dependence of galaxy bias as a function 
of luminosity and color. iTegmark et alj d2004ah applied a matrix- 
based method using pseudo-Karhunen-Loeve eigenmodes to mea- 
sure three power spectra from the SDSS galaxy distribution, allow- 
ing a quantification of the clustering anisotropy and a more accu- 
rate reconstruction of the real-space power spectrum than can be 
obtained from the angle-averaged redshift space power spectrum. 
Non-linear redshift space distortions, caused in part by the virial- 
ized motions of galaxies in their host dark matter haloes, create 
features known as Fingers-of-God (FOGs) alon g the line of sight 
in the redshift space galax y density field jDavis & Peebles!! 19831: 
iGramann et alj|l994l) . Both lTegmark et al.1 ( I2004ah and lCole et alj 
(2005) apply cluster-collapsing algorithms to mitigate the effects of 
FOGs before computing power spectra. Previous analyses have fit 
galaxy power spectra to linear {Percival et al. 2001 , 2007 ) or non- 
linear matter models dSpergel et alj2003l ; lTegmark et al.l2004l) . but 
did not attempt t o model the scale dependence of the galaxy bias. 
Cole ct al. I d2005h introduced a phenomenological model to account 
for both matter non-linearity and the non-trivial relation between 
the galaxy power spectrum P ga i(k) and matter power spectrum: 



P Ba i{k) 



Qk< 



1 + Ak 



Plin (&) • 



(l) 



where Pu n denotes the und erlying linear m atter power spectrum. 
For the 2dFGRS analysis, ICole et alj d2005l) fit A using mock 
galaxy catalogues and derive expected central values of Q. In the 
fit to the observed galaxy power spectrum, they allow Q to vary up 
to twice the expected value, which is supported by halo model cal- 
culations of the cosmological dependence of the galaxy P(k). This 
approach appears to work well for the case of 2dFGRS galaxies 
because it was calibrated on mock catalogues designed to match 
the properties of this galaxy popula t ion; h owever, its application 
to the LRG sample in lTegmark et alj |2006), where the best-fitting 
Q was much larger tha n for 2d FGRS galaxies, is questionable 
(seelReid et alj|2008l and lYoo et al]|2009l . but also lSanchez & Cold 
120081) . 

In this paper we focus our efforts on accurately modeling the 
relationship between the galaxy and matter density fields for the 
SDSS LRG sample. Several authors have studied this relation using 
the small and intermedia t e scale clustering in the SDSS LRG sam- 
ple dMasiedi et alj | 200tjl; IZehayi et alj|2005l: iKulkarni et alj|2007l : 
IWake et alj 120081; Izheng et alj 120081; Weid & Spergell l2009h and 
galaxy-galaxy lensing (Mandelbaum et al. 2006). The LRG selec- 
tion algorithm in the SDSS ( [Eisenstein et aU 2001) was designed to 
provide a homogenous galaxy sample probing a large volume with 
a number density, ulrg, which maximizes the effective survey 
volume V e ff(k) on the large scales of interest, k ~ .1 /t Mpc -1 . 
Vef f is given by dFeldman et alJll994l ; lTegmarklll997l) : 
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V Bf f{k) = 



n(r)P(fc) 



1 + n(r)P(k) 



(2) 



where P(k) denotes the measured galaxy power spectrum, n(r) 
the average galaxy number density in the sample at position r, 
and the integral is over the survey volume. The total error on 
P(k) is minimized (i.e., Veff is maximized) when nP ~ 1, 
which optimally balances cosmic variance and shot noise for a 
fixed number of galaxies. The LRG sample has proven its statisti- 
cal power through th e detection of the BAO ( Eise nstein et alj20o3 ; 
IPercival et al.l l2007t) . However, parameterizing the LRG power 
spectrum with a heuristic model for the non-linearity (Eqn.[TJ and 
marginalizing over fitting parameters limits our ability to extract 
the full cosmological information available from the power spec- 
trum shape and can introduce systematic biase s I Sanchez & Cold 



l2008t IDunklev et ai]|2009l : IVerde & . Peirisll2008l : iReid et al] 2008 ) . 

On sufficiently large scales, we expect galaxies to be lin- 
early biased with respect to the underlying m atter density field 
dMo & Whitell 19961 : IScherrer & Weinberg||l998l) . However, an of- 
ten overlooked consequence of a sample with tilrgPlrg ~ 1 
is that errors in the treatment of the shot noise can introduce sig- 
nificant changes in the measured shape of PLRo(k) and can be 
interpreted as a scale dependent galaxy bias. In the halo model pic- 
ture, the LRGs occupy massive dark matter haloes, which them- 
selves may not be Poisson tracers of the underlying matter density 
field, as they form a t the high peaks o f the initial Gaussian den- 
sity distribution (e.g. jBardeen et al.ll986h . Moreover, an additional 
shot noise-like term is ge nerated when multiple LRGs occupy indi- 
vidual dark matter haloes dPeacock & SnTitrteOOCtlCoorav & Shethl 
l2002h . Our approach is to first eliminate the one-halo contribution 
to the power spectrum by identifying groups of galaxies occupying 
the same dark matter halo, and then to calibrate the relation be- 
tween the power spectrum of the reconstructed halo density field, 
Phaio(k, p), and the underlying matter power sp ectrum, PpM(k) , 
using the A-body simulation results presented in lReid et alj ( 2008). 
As a result, the effects of non-linear redshift space distortions 
caused by pairs of galaxies occupying the same halo are dimin- 
ished. However, a further complication is that LRGs occupy the 
massive end of the halo mass function, and velocities of isolated 
LRGs within their host haloes could still be quite large. The details 
of the relation between LRGs and the underlying matter distribu- 
tion can then have a significant impact on the non-linear corrections 
to the power spectrum. 

The DR7 LRG sample has sufficient statistical power that 
the details of the relation between LRGs and the underlying mat- 
ter density field become important and need to be reliably mod- 
eled before attempting a cosmological interpretation of the data. 
This paper offers three sequential key improvements to the m odel- 
ing of LRG clusterin g compared with Eis enstein et al .1 d2005h and 
lTegmarketaT] < l2006h : 

• We reconstruct the underlying halo density fiel d traced by the 
LRGs before computing the power spectrum, while iTegmark et all 
(2006) apply an aggressive FOG compression algorithm. The re- 
constructed halo density field power spectrum deviates from the 
underlying matter power spectr um by < 4% at k = 0.2 h Mpc - , 
while the lTegmark et alj (120061) power spectrum differs by ~ 40% 
at k = 0.2 h Mpc" 1 dReid 



Mb ) power 
et alJ2 008) 



• We produce a large set of mock LRG catalogues drawn from 
A-body simulations of sufficient resolution to trace a halo mass 
range relevant to LRGs without significant errors in the small-scale 
halo clustering and velocity statistics (see Appendix A o flReid et alj 
2008). We present novel consistency checks between the mock and 



observed LRG density fields in halo-scale higher order clustering, 
FOG features, and the effective shot noise. 

• We use these tests along with the halo model framework to de- 
termine tight bounds on the remaining modeling uncertainties, and 
marginalize over these in our likelihood calculation. In contrast, 
lEisenstein et al. ( 2005) assume no uncertainty in their model LRG 
correlation function, an d lTegmark etafl (2006) marginalize over Q 
in Eqn.[T]with only an extremely weak prior on Q. 



This paper represents a first attempt to analyse a galaxy redshift sur- 
vey with a model that accounts for the non-linear galaxy bias and 
its uncertainty; other approaches that utilize the galax y distribu- 
tion r ather than the halo density field are in development dYoo et all 
120091) . 

In this paper we present and analyse a measurement of the 
power spectrum of the reconstructed halo density field from the 
SDSS DR7 LRG sample. DR7 represents a factor of ~ 2 increase 
in effe ctive volume over the anal yses presented in lEisenstein et al.1 
(2005) and lTegmark et alj (2006), and covers a coherent region of 
the sky. Section [2] describes the measurement of the reconstructed 
halo density field power spectrum, Phaio(k), along with the win- 
dow and covariance matrices used in our likelihood analysis. Sec- 
tion[3]describes the details of our model for the reconstructed halo 
power spectrum, Phaio(k, p). In Section|4]we summarize the tests 
we have performed for various systematics in our modeling of the 
relation between the galaxy and dark matter density field. We quan- 
tify the expected level of uncertainty through two nuisance param- 
eters and present several consistency checks between the model 
and observed reconstructed halo density field. In Section|5]we dis- 
cuss the cosmological constrain ts from Phaio(k) a lone as well as 
in combination wit h WMAP5 dPunklev et alj|2009l) and the Union 
supernova dataset dKowalski etalj|2008h . Section [ol compares our 
findings with the results of previous analyses of galaxy clustering, 
and Section|7]summarizes our conclusions. 

In a companion paper (iPercival et aljprert hereafter, P09) we 
measure and analyse BAO in the SDSS DR7 sample, of which 
the LRG sample considered here is a subset. BAO are detected 
in seven redshift shells, leading to a 2.7% distance measure at 
redshift 2 = 0.275, and a measurement of the gradient of the 
distance-redshift relation, this quantified by the distance ratio be- 
tween z — 0.35 and z — 0.2. We show in Section[5]that the results 
from these measurements are in agreement with our combined re- 
sults from BAO and the shape of the power spectrum calculated 
using just the LRGs. The results from these different analyses will 
be correlated because of the overlapping data used, so they should 
not be combined in cosmological analyses. The best data set to be 
used will depend on the cosmological model to be tested. While 
the inclusion of 2dFGRS and main SDSS galaxies in P09 provides 
a higher significance detection of the BAO, we show in Section 
15. 41 that the full power spectrum information provides tighter con- 
straints on both massive neutrinos and the number of relativistic 
species. 

Throughout the paper we make use of t wo specific cosmolog- 
ical models. The simulation set described in lReidetalJd2008l) and 
used to calibrate the model Phaio(k, p) adopts the WMAP5 recom- 
mended ACDM values: (Q. m , Q b , fi A , n s , <x 8 , h) = (0.2792, 0.0462, 
0.7208, 0.960, 0.817, 0.701). We refer to this model throughout the 
paper as our 'fiducial cosmological model.' To convert redshifts 
to distances in the computation of the Phaio(k), we adopt a flat 
ACDM cosmology with Q, m = 0.25 and fi A = 0.75. Throughout 
we refer to the power spectrum of several different density fields 



© 0000 RAS, MNRAS 000, 000-000 



4 Beth A. Reid et al. 



P(k) Definition Reference 



i'r,iu;(k) measured angle averaged redshift- space power spectrum of the LRGs 

Phalo(k) measured power spectrum of reconstructed halo density field 

Plin(k) linear power spectrum computed by CAMB Lewis et al. (2000) 

Pdm{^) theoretical real-space non-linear power spectrum of dark matter 

Pnw(k) theoretical linear power spectrum without BAO ("no wiggles") Eisenstein & Hu (1998) 

Pdamp{k) theoretical linear power spectrum with damped BAO (Eqn. llOt Eisenstein et al. (2007b) 

Phalo(k, p) model for the reconstructed halo power spectrum for cosmological parameters p Reid et al. (2008) 



Phalo,win(k, p) Phalo(k, p) convolved with survey window function (Eqn.|5J Percival et al. (2007) 

and directly compared with PhaloW m the likelihood calculation (Eqn.[6) 

Table 1. Definitions of power spectra referred to throughout the paper. 



and several theoretical spectra. Table Q] summarizes their defini- 
tions. 



2 DATA 

2.1 LRG sample 



The SDSS dYork et al.l l2000t) is the larg est galaxy survey 

etal.1 l200fj) 



ever produced; it used a 
to obtain imaging data in 



2.5m telescope l lGun 
5 passb ands u, g, 



i and z 



Fukugita et al. 1996; Gunn et al. 2006) . The images were reduced 
Stoughton etalJl2002l;|Pier et aU2003:llvezic et alj2004l) and cal- 



ibrated faogg et al I l200ll ; ISmith et alj|2002l ; iTucker et alj 12006) : 
IPadmanabhan et alj|2008l) . and galaxies were selected for follow- 
up spectroscopy. The second phase of the SDSS, known as SDSS- 
II, has recently finished, and the DR7 dAbazaiian etai]|2009T) sam- 
ple has recently been made public. The SDSS project is now contin- 
uing with SDSS-III where the extragalact ic component, the B aryon 
Oscillation Spectroscopic Survey (BOSS uSchlegel et alj2009l) , has 
a different galaxy targeting algorithm. DR7 therefore represents the 
final data set tha t will be released with the original target ing and 
galaxy selection (Eise nstein et alj200ll ; IStrauss et alj|2002h . 

In this paper we analyse a subsample containing 110 576 Lu- 
minous Red Galaxies (LRGs: lEisenstein et alj|200lh . which were 
selected from the SDSS imaging based on g, r and i colours, 
to give approximately 15 galaxies per square degree. The SDSS 
also targeted a magnitude limited samp le of galaxies for spec- 
troscopic follow-up dStrauss et al.l 120021) . The LRGs extend this 
main galaxy sample to 2 ~ 0.5, covering a greater volume. Our 
DR7 sample covers 7931 deg 2 (including a 7190 deg 2 contiguous 
region in the North Galactic Cap), with an effective volume of 
V e ff — 0.26 Gpc 3 h -3 , calculated with a model power spectrum 
amplitude of 10 4 /i _3 Mpc 3 . This power spectrum amplitude is ap- 
proximately correct for the LRGs at k ~ 0.15 h Mpc~ . For com- 
pariso n, the effective volume of the sample used bv lEisenstein et al] 
(l2005l) was V„ff = 0. 13Gpc 3 h- 3 , and V eff = 0.16Gpc 3 h" 3 



Tegma rk et al. I d2006h : this work represents a factor of ~2 in- 



crease in sample size over these analyses. The sample is the same as 
that us ed in P09, and its construction follows that of IPercival et al] 
d2007f) . albeit with a few improvements. 

We use SDSS Galactic extinction-corrected Petrosian 
magnitudes calibra t ed u sing the "ubercalibration" method 
(Padmanabhan et al. 2008). However, we find that the power 
spectrum does not change signi ficantly when on e adopts the old 
standard calibration instead dTuckeret afl l2006t). Luminosities 
are K-corrected using the methodology of Blanton e t al] d2003h : 
iBlanton et al] ]2003bl) . We remove LRGs that are not intrinsically 
luminous by applying a cut Mo.i r < —21.8, where Mo.i r is our 
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Figure 1. Fits to the redshift distributi ons for the LRG selection used in this 
work (solid cur ves) and fhelZehavi et all j2005t) -23.2 < M g < -21.2 
sample used in lTegmark et all (2006) (dashed curves). Upper panel: n(z) 
vs redshift in units of 10~ 4 (h^ 1 Mpc)~ 3 Lower panel: N(< z) 
1' dzn(z)dV/dz (arbitrary overall normalization). 



estimate of the absolute magnitude in the r-band for a galaxy at 
z = 0.1. 

Spectrosco pic LRG targets were selected using two color- 
magnitude cuts dEisenstein et al] 12001). The tiling al gorithm en- 
sures nearly complete samples telanton et all [2003). However, 
spectroscopic fiber collisions prohibit simultaneous spectroscopy 
for objects separa ted by < 55", leavin g ~ 7% of targeted objects 
with out redshifts dMasiedi et al.| [2006). We correct for this effect 
as in lPercival et ali i 20071) : foran LRG lacking a spectrum but 55" 



from an LRG with a redshift, we assign both galaxies the measured 
redshift. If the LRG lacking a redshift neighbors only a galaxy from 
the low redshift SDSS main sample, we do not assign it a redshift. 
These galaxies are assumed to be randomly distributed, and simply 
contribute to the analysis by altering the completeness, the fraction 
of targeted galaxies with good redshifts, in a particular region. The 
impact of the fiber collision correction is addressed in Appendix 
[B3]and Appendix lB4l 

Fig. Q] compares the number density as a function of red shift 
for the LRG s election used in this p aper dPercival et al]|2007l) and 
the on e used in[Tegmark et ai]( l2006l) and presented in lzehavi et all 
d2005l) . The main differences are that our selection includes a small 
number of galaxies at z < 0.15, and our cut on the intrinsic lumi- 
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nosity of the LRGs slightly reduces the number density of galax- 
ies at high 2. The different selections produce a similar number of 
galaxies per unit volume, and we expect no difference between the 
samples on the large scale structure statistics of interest here. 

2.2 Recovering the halo density field 

In real space, the impact of more than one LRG per halo on the large 
scale power spe ctrum can be accurately modeled as an additional 
shot noise term dCoorav & Shethl2002l : lReid et alj2008l) . However, 
this picture is much more complicated in redshift-space because 
of the velocity dispersion o f the LRGs shift s them along the line 
of sight by ~ 9 h' 1 Mpc l iReid et alj|2008h . and the distribution 
of intrahalo velocities has long tails. This shifting causes power 
to be shuffled between scales and causes even the largest scale 
modes along the line of sight to be damped by these FOG fe atures 
jDavis & Peebleslll983l ; |Peacock & Doddsll 1 9941 : 1 S eli akll200 lb . We 
substantially reduce the impact of these effects by using the power 
spectrum of the reconstructed halo density field. 

We follow the Counts-In-Cylinders (CiC) technique in 
Rcid et al. | d2008h to identify LRGs occupying the same halo and 
thereby estimate the halo density field. Two galaxies are consid- 
ered neighbors when their transverse comoving separation satisfies 
Ar± < 0.8 hT 1 Mpc and their redshifts satisfy Az/(1 + z) < 
Av p /c = 0.006 (Av p = 1800 km s -1 ). A cylinder should be a 
good approximation to the density contours of satellites surround- 
ing central galaxies in redshift space, as long as the satellite ve- 
locity is uncorrelated with its distance from the halo centre and 
the relative velocity dominates the separation of central and satel- 
lite objects along the line of sight. Galaxies are then grouped with 
their neighbors by a Friends-of-Friends (FoF) algorithm. The re- 
constructed halo density field is defined by the superposition of the 
centres of mass of the CiC groups. We refer to the power spectrum 
of the reconstructed halo density field as Phaio(k); it is our best es- 
timate of the power spectrum of the haloes traced by the LRGs. For 
comparison we also compute the power spectrum without applying 
any cluster-collapsing algorithm, PLRo(k). 

Our reconstructed halo density field contains 104 337 haloes 
derived from 110 576 LRGs. 



2.3 Calculating power spectra, window functions and 
covariances 

In this paper we focus on using the angle-averaged power spectrum 
to derive constraints on the underlying linear theory power spec- 
trum. On linear scales the redshift space power spectrum is pro- 
portio nal to the real space power spectrum dKaise J 1 9871 : iHamiltonl 
1998). Our halo density field reconstruction mitigates the effects of 
FOGs from objects occupying the same halo. Though we do not 
explore it here, we expect that our halo density field reconstruction 
will be useful to an analysis of redshift-space anisotropics (e.g., 
lHatton&Colelll999h . 

The methodology for calculating the power spectrum of the 
reconstructed ha lo density field , Phaio(k), is based on the Fourier 
method of iFeldman et al.l dl994l) . The halo density is calculated by 
throwing away all but the brightest galaxy where we have located 
a set of galaxies within a single halo. This field is converted to 
an over-density field by placing the haloes on a grid and subtract- 
ing an unclustered "random catalogue", which matches the halo 
selection. To calculate this random catalogue, we fit the redshift 
distributions of the halo sample with a spline model dPress et all 




0.05 0.1 0.15 0.2 

k / h Mpc" 1 



Figure 2. Top panel: Measured Phalo(k) bandpowers. Error bars indicate 
%/ Ca (Eq. [5J. Middle panel: Correlations between data values calculated 
using Log-normal catalogues, assuming our fiducial cosmological model. 
Bottom panel: The normalized window function for each of our binned 
power spectrum values with 0.02 < k < 0.2 h Mpc -1 . Each curve shows 
the relative contribution from the underlying power spectrum as a function 
of k to the measured power spectrum data. The normalisation is such that 
the area under each curve is unity. For clarity we only plot curves for every 
other band power. 



fl992h (shown in Fig.Q3, and the ang ular mask was dete rmined us- 
ing a routine based on a HEALPIX ( Gorski et al. 2005) equal-area 
pixelization of the sphere as in dPercival et al.ll2007l) . This proce- 
dure allows for the variation in radial selection seen at z > 0.38, 
which is caused by the spectroscopic features of the LRGs moving 
across the wavebands used in the target selection. The haloes and 
randoms are weighted using a luminosity-dependent bias model 
that normalizes the fluctuations to the amplitude of L* galaxies 
dPercival et ai]|2004h . To do this we assume that each galaxy used 
to locate a halo is biased with a linear deter ministic bias model, and 
tha t this bias depends on Mo.i r according to lTegmark et al.l d2004a) 
an d lZehavietal1 d2005), where Mo.i r is the Galactic extinction and 
K-corrected r-band absolute galaxy magnitude. This procedure is 
similar to that adopted by P09. 

The power spectrum was calculated using a 1024 3 grid in a 
series of cubic boxes. A box of length 4000 Mpc was used 
initially, but we then sequentially divide the box length in half and 
apply periodic boundary conditions to map galaxies that lie out- 
side the box. For each box and power spectrum calculation, we 
include modes that lie between 1/4 an d 1/2 the Nyquis t frequency 
(similar to the method described by ICole et all 120051) . and cor- 
rect for the smoothing effect of the cloud-in-cell assignment used 
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to locate galaxies on the grid (e.g. iHocknev & Eastwood I [l98ll 
chap. 5). The power spectrum is then spherically averaged, leav- 
ing an estimate of the "redshift-space" power. The upper panel of 
Fig. [2] shows the shot-noise subtracted bandpowers measured from 
the halo density field, calculated in bands linearly separated by 
Afc = 0.004 /iMpc" 1 . This spacing is sufficient to retain all of 
the cosmological information. 

The calculation of the likelihood for a cosmological model 
given the measured bandpowers Phaloik) requires three additional 
components determined by the survey geometry and the proper- 
ties of the galaxy sample: the covariance matrix of measured band- 
powers dj, the window function W (ki,k n ), and the model power 
spectrum as a function of the underlying cosmological parameters, 
Phaio(k, p). The calculation of model power spectra is considered 
in Section[3] 

The covariance matrix and corresponding correlation coeffi- 
cients between bandpowers i and j are defined as 



corr(ki,kj) — 



Cij = {Phalo{ki)Phalo{kj)) - {Phalo{ki)){Phalo{kQ) 

(4) 



The covaria nce matrix was calcula t ed from 10 4 Log- Normal (LN) 
catalogues dColes & Jonesl Il99ll ; ICole et alj l2005h . Catalogues 
were calculated on a (512) 3 grid with box length 4000 h' 1 Mpc 
as in P09, where LN catalogues were similarly used to estimate 
covariance matrices. Unlike iV-body simulations, these mock cat- 
alogues do not model the growth of structure, but instead return 
a density field with a log-normal distribution, similar to that seen 
in the real data. The window functions for these catalogues were 
matched to that of the halo catalogue. The input power spectrum 
was a cubic spline fit matched to the data power spectra, multi- 
plied by a damped ACDM BAO model calculated using CAMB 
ilLewis et alhoOOl) . The recovered LN power spectra were clipped 
at 5a to remove extreme outliers which contribute less than 0.05% 
of the simulated power spectra, and are clearly non-Gaussian. This 
covariance matrix calculation matches the procedure adopted by 
P09. The middle panel of Fig. [2] shows the correlations expected 
between band-power s calculated using this procedure. 

As described in ICole et al.l J2005t) . the window function can 
be expressed as a matrix relating the theory power spectrum for 
cosmological parameters p and evaluated at wavenumbers k n , 
Phaio(k„,p), to the central wavenumbers of the observed band- 
powers kc 

Phalo,win(ki,p) = ^2 W(ki, k„)P ha lo(kn,p) ~ W(&i, 0). (5) 

71 

The term W(ki, 0) arises because we estimate the average halo 
density from the sample, and is related to the integral constraint in 
the correlation function JPercival et all [2007). The window func- 
tion allows for the mode-coupling induced by the survey geometry. 
Window functions for the measured power spectrum (Eqn. 15 of 
IPercival eta l, 2004) were cal culated as described in lpercival et al.l 
d200ll) . ICole et aljj2005h . and lPercival et al.ld2007l) : an unclustered 
random catalogue with the same selection function as that of the 
haloes was Fourier transformed using the same procedure adopted 
for our halo overdensity field described above. The shot noise was 
subtracted, and the power spectrum for this catalogue was spheri- 
cally averaged, and then fitted with a cubic spline, giving a model 
for W (ki, k n ). For ease of use this is translated into a matrix by 
splitting input and output power spectra into band powers as in 
Eqn. |U 



The window functions W(ki, kj) and the corresponding cor- 
relation coefficients for every other bandpower are shown in the 
lower panel of Fig. [2] In addition to the window coupling for 
nearby wavenumbers, there i s a beat-couplin g to survey-scale 
modes ^Hamilton etal] 12004 feeid et al.l l2008h . That is, density 
fluctuations on the scale of the survey couple to the modes we 
can measure from the survey. However, this effect predominantly 
changes only the amplitude of Phaloik), which is marginalized 
over through the bias para meter bn m Eqn. | 15| below. Fig.|2]can be 
compared with Fig. 10 in lPercival et al.l J2007l) . where the windows 
and correlations were presented for the SDSS DR5 data. For the 
DR5 plot, variations in the amplitude were removed leaving only 
the small-fc difference couplings. The power spectrum, window 
functions, and inverse covariance matrix are electronically avail- 
able with the likelihood code we publicly release (see Section[5). 



2.4 Phaioik) likelihood 

We assume that the likelihood distribution of the power spectrum 
band powers is close to a standard multi-variate Gaussian; by the 
central limit theorem, this should be a good approximation in the 
limit of many modes per band. The final expression for the likeli- 
hood for cosmology p is then 



2 In L(p)= X 2 (p) = ^A ! C7r 1 A J 



(6) 



where Aj = [(Phaio(ki) - Phaio,win(ki, p)] . 

A single comoving distance-redshift relation Xfidiz), that of 
a flat, Q m — 0.25 cosmology, is assumed to assign positions to the 
galaxies in our sample before computing Phaio(k). Rather than re- 
computing Phaio(k) for each com oving distance-redshift relation 
to be tested, IPercival et ail (2007) and P09 account for this when 

nodels by 
~^]|2005h ( 

tifies the model dependence of the conversion between (ra, dec, z) 
and comoving spatial coordinates when galaxy pairs are distributed 
isotropically: 

r 1 V3 

(l + z) 2 D A (z) 



evaluating the likelihood of other c osmological models by alter- 
ing the window function. Dy(z, p) dEisenstein et alj [2005) quan- 



D v (z) 



H[z) 



(7) 



wh ere Da(z) is the phy sical angular diameter distance. Follow- 
ing iTegmark et al.l J2006h we partially correct for the discrepancy 
between the fiducial model Xfid(z) and the x( z ) °f me model to 
be tested by introducing a single dilation of scale. To first order, 
changes in the cosmological distance-redshift model alter the scale 
of the measured power spectrum through Dv(z), so we introduce 
a scale parameter that depends on this quantity, 



a s d{z) 



Dv{z) 

£)fiducial( 2 ) ■ 



(8) 



Strictly, we should allow for variations in a sc i across the redshift 
range of the survey, as in P09. However, to first approximation we 
can simply allow for a single scale change at an effective redshift 
for the survey Z e ff- When comparing Phaio(k), computed using 
Xfid(z), with a model comoving distance-redshift relation x{ z , p), 
in practice we use0 

1 This correction was incorrectly applied in previous versions of COS- 
MOMC, and is corrected in the code we release. This correction is primarly 
important for constraining the BAO scale rather than the turnover scale, and 
so previous analyses with COSMOMC should be minimally affected. 
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Aj = [Phalo(ki 



Phc 



In Appendix I A2I we verify that this approximation is valid for our 
sample with z e ff = 0.313. 

In our cosmological analysis, we include modes up to k max = 
0.2 ftMpc -1 , where the model power spectrum deviates from the 
input linear power spectrum by < 15%. We also impose a conser- 
vative lower bound at k m i n — 0.02, above which galacti c extinc- 
tion corrections (see the analysis in IPercival et alj|2007h . galaxy 
number density modeling, and window function errors should be 
negligible. 

P09 present a detailed analysis demonstrating that the BAO 
contribution to the likelihood surface is non-Gaussian; this is in 
large part due to the relatively low signal-to-noise ratio of the BAO 
signature in our sample. Therefore, to match expected and recov- 
ered confidence intervals, P09 find that the covariance matrix of 
the LRG-only sample must be inflated by a factor 1.1 2 = 1.21. 
Though our likelihood surface incorporates constraints from the 
shape of the power spectrum, for which the original covariance 
matrix should be accurate, we conservatively multiply the entire 
covariance matrix by this factor required for the BAO constraints 
throughout the analysis. Therefore our constraints likely slightly 
underestimate the true constraints available from the data. This fac- 
tor is already included in the electronic version we release with the 
full likelihood code. 



3 MODELING THE HALO POWER SPECTRUM 

We consider three effects that cause the shape of Phalo(k, p) to 
deviate from the linear power spectrum, Pu n (k, p), for cosmolog- 
ical parameters p. We will assume that these modifications of the 
linear power spectrum can be treated independently. These effects 
are the damping of the BAO, the change in the broad shape of the 
power spectrum because of non-linear structure formation, and the 
bias because we observe galaxies in haloes in redshift space rather 
than the real space matter distribution. We also need to consider the 
evolu tion of the s e effe cts with redshift. 

iReid et al.l J2008I) construct a large set of mock LRG cata- 
logues based on A^-body simulations evaluated at a single cosmo- 
logical model Pfid- We use these catalogues to calibrate the model 
halo power spectrum, and make detailed comparisons between the 
observed and mock density fields in Appendix |B1 



3.1 BAO damping 

The primary effect of non-linear structure format ion and peculiar 
velociti es on the BAOs is to damp them at large fc. lEisenstein et al.l 
(2007b) showed that this can be accurately modelled as a Gaussian 
smoothing, where 



(9) or a tracer like the LRGs. We fix Ohaio, i.e., the value of a ap- 
propriate for the reconstructed halo density field, using fits to the 
reconstructed halo den sity field power s pectrum in the mock LRG 
catalogues presented in lReid etal and shown here in Fig. [4] 

We performed tests in the ACDM case which demonstrate that cos- 
mological constraints are not altered when a is allo wed to vary with 
cosmol ogy p according to the dependence given in lEisenstein et al.l 
J2007bh . and in Appendix lA3l we show that using a spline fit to Pu n 
instead of the Eisenstein & Hu ( 1998) formula for P nw does not 
affect the likelihood surface in the region of interest. 



3.2 Non-linear structure growth 

As the small perturbations in the early universe evolve, gravita- 
tional instability drives the density field non-linear, a nd power on 
small scales is enhanced as structures form. HALOFIT dSmith et"al] 
120031) provides an analytic formalism to estimate the real space 
non-linear matter power as a function of the underlying linear mat- 
ter power spectrum. While Eqn.[l0]accounts for the effects of non- 
linear growth of structure on the BAO features in Phalo(k, p), 
HALOFIT provides a more accurate fit to the smooth component 
of the non-linear growth in the quasi-linear regime (k ^ 0.2) when 
evaluated with an input spectrum P nw (k, p) rather than the linear 
matter power spectrum containing BAO wiggles: 



damp 



(fc,p,a) = Pi in (fc,p)e 2 +P n 



r(*,p) (l 



Here Pii n(fc, p) is the line ar matter power spectrum computed by 
CAMB dLewis et alj feoOO) and shown in the upper left panel of 
Fig. [3] for our fiducial cosmolo gical model. P nw (fc,p), defined 
by Eqn. 29 of lEisenstein & Hul dl998l) . is a smooth version of 
Piin(fc, p) with the baryon oscillations removed. The upper right 
panel of Fig. [3] shows the ratio Pu n (k) / P n w(k) for our fiducial 
cosmology. The amplitude of the damping is set by a and depends 
on the cosmological parameters, whether the power spectrum is in 
real or redshift space, and whether we are considering the matter 



rhaiofit(k,p) = 



halofrt.nw 



(fc,p) 



Pnw(k, p) 

PDM,halofit{k,p) = Pdamp{k,p,a m )rhalofit(k,p) 



(ID 

(12) 



Eqn. [12] is our modified HALOFIT model real space power spec- 
trum, using Eqn.QJJ]to account for BAO damping and HALOFIT for 
the smooth component. The bottom left panel of Fig. [3] shows that 
PDM{k)/P damp {k, <7 m ) and TWofit agree at the ~ 1.5% level for 
k ^ 0.2 in our fiducial cosmology. Since we normalize the final 
model Phaio(k, p) using our mock catalogues at the fiducial cos- 
mology Pfid, in practice HALOFIT only provides the cosmological 
dependence of the non-linear correction to the matter power spec- 
trum: 



rDM,damp(k,p) — 



rhaiofit(k, p) 



r(k, pad) 



rhalofit(k, Pfi d ) Pdamp(fc, Pfid, ODm) 



(13) 



•(10) 



rDM,dam P (k, p) is our model for the ratio of the non-linear matter 
power spectrum to the damped linear power spectrum. The nor- 
malization of roM,dam P accounts for the small offset between the 
A-body and HALOFIT results in Fig. [3] at the fiducial cosmology. 
In the space of cosmologies consistent with the data, the small 
cosmology-dependence of this correction is primarily through a$. 
In Section \52\ we find that the LRG-only likelihood surface is in- 
dependent of the assumed value of a% over the range 0.7 to 0.9. 



3.3 Halo bias 

In our likelihood calculation we marginalize over the overall am- 
plitude of Phaio(k), so in this Section we are concerned only with 
the scale dependence of th e relation between th e reconstructed halo 
and matter power spectra. ISmith et alj ( 120071) show that the scale 
dependence of halo bias in real space is large for the most mas- 
sive haloes, but should be rather wea k for the halo ma ss range 
which host the majority of the LRGs; Matsubara] d2008f) demon- 
strates t his analytically in redshift space in the quasi-linear regime. 
Indeed, [Reid et alj d2008l) find that the power spectrum of the (red- 
shift space) reconstructed halo density field is nearly linearly biased 
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Figure 3. Upper left panel: Power spectra for the fiducial cosmology. 
The solid curve is Pu n (k) and the dashed curve is Pnw(fc)fhalofit> 
the non-linear power spectrum from HALOFIT using -Pnw 

(fc) as 

the input. Upper right panel: F\ ln (k) / P n w (fc) ■ Bottom left panel: 
PpM(k) I Pdamv (k, o~m) measure d in TV-body simulation snapshots at 
zmid, reported in Reid etal. (2008), compared with the smooth correction 
predicted by HALOFIT, r haloflt . Bottom right panel: r haloflt at {z NEAR , 
zmid, z FAR } = {0.235, 0.342, 0.421}. 



with respect to the underlying real space matter power spectrum for 
fc < 0.2 h Mpc" 1 and our fiducial ACDM model, and we assume 
this should remain approximately true in the narrow range of cos- 
mologies consistent with the data. For the fiducial cosmology, we 
can use our simulations to calibrate the relation between the halo 
and matter spectra: 

Pfid)/ Pdamp{k, p, (Thalo) 



rhalo,Du{h, Pfid) 



(14) 



Pfid) / ' Pdamp{k, p, (7m) 

This is our model for the smooth component of the bias between the 
halo and dark matter power spectra. To account for any dependence 
of rhaio.DM(k, Pfid) on the cosmological model and other remain- 
ing modeling uncertainties, we introduce a smooth multiplicative 
correction to the final model Phalo(k, p) containing three nuisance 
parameters bo, ai and a?,: 

F nuis (k) = b 2 + ai (A) + a2 (j-)'^ (15) 

where we set fc* = 0.2 h Mpc -1 . The parameter bo is the effective 
bias of the LRGs at the effective sample red shift, z e ff, relative to 
L* galaxies (Eqn. 18 of IPercival et al.ll2004) . In Section|4]we will 
use consistency checks between the observed and mock catalogue 
galaxy density fields as well as the halo model framework to es- 
tablish the allowed region of ai — ai parameter space. An allowed 
trapezoidal region in ai — ai space is completely specified through 
two parameters, Uo.i and no. 2- These two parameters specify the 
maximum absolute deviation allowed by F nU i S (k)/bo away from 
1 for fc ^ 0.1 (uo.l) and 0.1 ^ fc ^ 0.2 (110.2). When evaluating 
the likelihood of a particular cosmological model we marginalize 
analytically over bo using a flat prior on bo 0, and we marginalize 
numerically over the allowed ai — ai region with a flat prior in this 
region. We discuss the impact of these priors on the cosmological 
constraints in Appendix [C] 



3.4 Model fits and evolution with redshift 

Our final model halo power spectrum at fixed redshift treats each 
of the three non-linear effects independently: Eqn. QJJ] converts 
the linear power spectrum to the damped linear power spectrum, 
r dm, damp converts the damped linear power spectrum to the real 
space non-linear matter power spectrum, r na io.DM converts the 
real space non-linear matter power spectrum to the redshift space 
reconstructed halo density field power spectrum (assuming this re- 
lation is cosmology independent), and F nU i 3 (fc) allows for smooth 
deviations from our model due to modeling errors, uncertainties, 
and unaccounted cosmological parameter dependencies: 



Phalo(k,p) = Pdam P (k,p)r D M,da m p(k,p) X 
Thalo,DM{k, Pfid)Fnuis(k). 



(16) 



For this multiplicative model, the PD!vi(fc, Pfid )/P 'da m P (k, p, a m ) 
terms from Eqns.[T3]and fj4]cancel, so calibration of the model only 
requires fits to er halo and Ph a io(k, Pfid)/Pdam P (k, Pfid, o^alo) 
using the mock catalogues. 

The model in Eqn. [16] is strictly only valid at a single red- 
shift. In order to match our model to the observed redshift dis- 
tribution of the LRGs and their ass ociated haloes, we use the 
mock halo catalogues constructed in iReid et alj d2008l) at three 
redshift snapshots. These are centered on the NEAR (zmear = 
0.235), MID (zmi d = 0.342) and FA R (z FAR = 0.421) 
LRG subsamples of iTegmark et alj J2006T) . Fig. [4] shows our fits 
to Phaio(k,p f id)/PH n {k,p S id) for each redshift snapshot. We 
first fit for <Jhaio in Eqn. 1 101 using our LRG mock catalogue re- 
sults Phaio{k, Pfid)- We include modes between k — Oh Mpc -1 
and fc = 0.2/i Mpc -1 in the fit and marginalize over an arbi- 
trary fourth order polynomial to account for the smooth devia- 
tions from Pdamp with fc. We find a ha to, near = 9.3 /i" 1 Mpc, 
(Jhaio, mid = 9.2 h' 1 Mpc, and a hato , FAR = 9.2 h' 1 Mpc. 
T hese numbers are rough ly consistent with the results presented 
in lEisenstein et al.1 J2007bl) . and are somewhat degenerate with the 
smooth polynomial correction. 

After fixing these values for <Jhaio, we calibrate the smooth 
component of the model, rDM,da m p(k,Pfid)r ha io,DM(k,pfid). 
For fc ^ 0.2 we fit P ha io(k, p/ id )/Pdamp(fc, Pfid, &haio) to a sec- 
ond order polynomial, and a fourth order polynomial for fc 5j 0.5. 
This component of the fit is shown in the first three panels of 
Fig. E] by the dotted curves, while the solid lines show the full 
fit to P ha io(k, p fid) /Pu n (k, pfid). Both the BAO-damping and 
smooth increase in power with fc are well described by our fits out 
tofc = 0.5/i Mpc" 1 . 

Our final model for the reconstructed halo power spectrum is 
a weighted sum over our model Phaio(k, p) (Eqn. 1 1 6b from each 
of the NEAR, MID, and FAR redshift slices fit in Fig.|4] 



Phalo(k,p) = 



-NEAR,AIID,FAR 



WiPhalo(k,p, Zi) , 



(17) 



where Wi specifies the weight of each redshift subsample. The 
lower right panel of Fig.[3]shows that the smooth correction for the 
non-linear matter power spectrum varies by < 1% over the redshift 
range of the LRGs. Moreover, the lower right panel of Fig.|4]shows 
that the relative shape of the power spectrum of the reconstructed 
halo density field varies by ± ~ 2.5% between the redshift sub- 
samples, so moderate biases in the determination of these weights 
will induce negligible changes in the predicted shape Phaio(k, p). 

In the limit that most pairs of galaxies contributing power to 
mode k come from the same redshift, the fractional contribution to 
the power spectrum from a large redshift subsample is 
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Figure 4. BAO-damping times polynomial fits to 
Phalo(k,Pfid)/Plin(k,Pfid) for o™ mock NEA R, MID, and FAR 
LRG reconstructed halo density field subsamples in Reid et al] j2008l) ; 
{znear, zmid, zfar} = {0.235, 0.342, 0.421}. The smooth compo- 
nent of these fits (dashed curves) enter our model Phalo(k,p) through 
Eqns. d!4t . while the amplitude of the BAO suppression &\ alo enters in 
Eqn. llQI Lower right panel: Ratio of the shape of the smooth components 
for the NEAR and FAR redshift subsamples to the MID redshift subsample. 



w(Zmin, Zmax) / U (z) ——(Iz , (18) 

J Zmin b 2 {z) dz 

where n(z), b(z), and w(z) specify the average number den- 
sity, bias, and weig ht of the sample at redshift z as defined in 
IPercival et al.l d2004l) . Since the integrand is slowly varying with 
redshift, this approximation should be fairly accurate. We derive 
weights wnear = 0.395, wmid = 0.355, and wfar = 0.250. 

3.5 Comparison with fiducial model Phaio(k, Pad) 

Our fiducial Phaio(k,p) model is calibrat ed on simulations wit h 
the WMAP5 recommended parameters dKomatsu et all 120091) : 
(tt m ,Q b ,n A ,n s ,a a ,h) = (0.2792, 0.0462, 0.7208, 0.960, 0.817, 
0.701). For the 45 observed bandpowers satisfying 0.02 < k < 
0.2/iMpc -1 , x 2 ~ 44.0 if we hold nuisance parameters a± = 
a,2 = and choose bo to minimize x' 2 ', our fiducial model is there- 
fore sufficiently close to the measured Phaio(k) to be used to cal- 
ibrate the cosmology-dependent model. The best-fitting nuisance 
parameters within the allowed range that we determine in Section 
1431 en = 0.172 and a 2 = -0.198, lower the X 2 to 40.9 for 42 
DOR The best-fitting model to the LRG-only likelihood presented 
in Section |5~T| is lower by only Ax 2 ~ 1.7 for the same treatment 
of the three nuisance parameters. 



4 QUANTIFYING MODEL UNCERTAINTIES AND 
CHECKS FOR SYSTEMATICS 

While the non-linear evolution of a collisionless dark matter den- 
sity field can be accurately studied using A^-body simulations, there 
remain many uncertainties in the mapping between the galaxy and 



matter density field. We first review the generic halo model predic- 
tions for a galaxy power spectrum, which provide the context for 
exploring the uncertainties in the relation between the galaxy and 
matter density fields. We summarize the results of Appendix [B] 
which presents our modeling assumptions and consistency checks 
between the mock catalogue and SDSS DR7 LRG density fields 
that constrain the level of deviation from our modeling assump- 
tions. The ultimate goal of this Section is to establish physically- 
motivated constraints on the nuisance parameters a\ and a 2 in 
Eqn.[l5]by determining urj.i and no. 2 defined in Section [331 These 
nuisance parameter constraints will then be used to compute cos- 
mological parameter constraints in Section|5] 

4.1 Galaxy power spectra in the halo model 

In the simplest picture for a galaxy power spectrum in the halo 
model, one considers a separation of the pairs into galaxies occu- 
pying the same dark matter halo, which contribute to P lh (k), and 
those oc cupying different dark matter haloes, which contribute to 
P 2h (k) JCooray & Shethl|2002h : 

P 3 al(k) = PlUk) + Pgal(k) (19) 

f# = [dMn(M) {N ^ N ^- 1)lM)) (20) 

J n gal 

P 2 9 Uk) = b 2 gal P DM (k). (21) 

On large scales, treating the haloes as linear tracers of the under- 
lying matter density field (Eqn. l21b and ignoring the spatial exten t 
of haloes in Eqn. |20] are good approximations dReidetal.H 2008). 
Therefore, in real space, the dominant effect of the inclusion of 
satellite galaxies is an excess shot noise given by Eqn. |20] though 
they also upweight highly biased halo pairs and slightly increase 
bgai as well. However, in redshift space, satellite galaxies are sig- 
nificantly displaced along the line of sight from their host haloes 
by the FOGs, and power is shuffled between scales, and even the 
largest scale modes along the line of sight are damped by the FOG 
smearing. There will be residual non-linear redshift space distor- 
tions in the reconstructed halo density field from imperfect recon- 
struction, and potentially from peculiar motion of isolated LRGs in 
their host haloes as well. 

4.2 Summary of tests for systematics and remaining 
uncertainties 

In the context of the halo model, both uncertainty in the distribu- 
tion of galaxies in groups as it enters Eqn.|20]and uncertainty in the 
structure of the FOG features will introduce uncertainty in the rela- 
tion between the reconstructed halo and matter density fields, and 
thus their power spectra. Appendices IB ll and lB2l d iscu s s the m odel- 
ing assumptions we have used to derive the lReid et al. mock 
LRG catalogues from A^-body simulation halo catalogues, and state 
the expected impact on the relation between the reconstructed halo 
and matter power spectra. 

Appendix lB 3 l introduces several distinct consistency checks of 
the uncertainties in Appendices IB 1 1 and IB2I In Section l2~2l we de- 
fine the CiC group finder by which we identify haloes. We demon- 
strate that this group finder produces group multiplicity functions 
that are in good agreement between the mock and observed LRG 
density fields, once fiber collisions are accounted for. While this 
agreement demonstrates that our mock catalogues reproduce small 
scale higher-order clustering statistics and FOG features of the ob- 
served density field, this is not a consistency check since the mocks 
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were designed to match these statistics. We find consistency when 
we compute a second CiC group multiplicity function allowing a 
wider separation between pairs perpendicular to the line of sight 
(Ar± = 1.2 ft" 1 Mpc). If the observed satellite galaxies were sig- 
nificantly less concentrated than in our mock catalogues, we would 
detect these galaxies when Ar± increases from 0.8 ft" 1 Mpc to 

1.2 ft" 1 Mpc. From this comparison we conclude that residual shot 
noise errors from inaccurate halo density field reconstruction are 
~ 2% of the total shot noise correction and do not dominate our 
systematic uncertainty. The second consistency check between the 
mock and observed LRG catalogues is the distribution of line of 
sight separations between pairs of galaxies in the same CiC group 
(Fig |B2t . This check probes the accuracy of our model of the FOG 
features coming from galaxies occupying the same halo, and the 
agreement we find indicates that the residual FOG features in the 
reconstructed observed and mock halo density fields will be in 
satisfactory agreement. Appendix IB4I presents the difference be- 
tween the power spectra with and without the halo density field 
reconstruction preprocessing step (Phaio{k) and Plug (k), respec- 
tively). This difference agrees with the mock catalogues, provided 
one carefully accounts for the impact of fiber collisions. In other 
words, while the treatment of fiber collisions can substantially im- 
pact PLRc(k), Phaio(k) is unaffected. In Appendix lB5l we demon- 
strate that the luminosity weighting used to compute Phaio(k) but 
not accounted for in the mock catalogues does not alter the effec- 
tive shot noise level of Phaio(k). Appendix |B6] presents evidence 
that the cosmology dependence of the model Phaio(k,x>) i s suf- 
ficiently accu rate. Finally, we n ote that lLunnan e t al. (prep) have 
compared the lReid et alj J2008h mock catal ogue genus curve with 
the observed genus curves dGott e t al. 2009), and find good agree- 
ment with no free parameters. 

As discussed in detail in Appendix IB2I the vast majority of 
LRGs (~ 94%) ar e expected to reside at the centre o f their host 
dark matter haloes feheng et al.ll2008l : iReid et alj|2008h . The prin- 
cipal modeling uncertainty we identify in Appendix|B]is the veloc- 
ity of these central LRGs within their host haloes; substantial in- 
trahalo velocities for these galaxies will suppress power in a scale- 
dependent manner (Fig. lBU . Note that none of the tests from Ap- 
pendix [B] can directly constrain the level of central LRG velocity 
dispersion. 

4.3 Constraints on F nu i a (fe) 

In Section [3731 we introduced a quadratic function F nu is{k) to 
account both for errors in our modeling at the fiducial cosmol- 
ogy and for any errors in the cosmology dependence of our 
model. We parametrized the amplitude of the total modeling un- 
certainty through uo.i and tto.2- These parameters, which we de- 
termine in this subsection, specify the maximum fractional devi- 
ation from the model power spectrum at k — 0.1 ft Mpc" 1 and 
k — 0.2 ft Mpc -1 , respectively. We choose these values of k be- 
cause k ^ 0.1 is usually considered safely in the linear regime, 
while k = 0.2 ft Mpc -1 is the maximum wavenumber we attempt 
to model. 

The dominant uncertainty in our model is in the relation be- 
tween the power spectrum of the reconstructed halo density field 
and the underlying matter power spectrum, which we describe by 
Eqn. 1141 At k = 0.1 ft Mpc -1 in the mock catalogues, the re- 
constructed halo density field and the redshift space central galaxy 
power spectra agree well below the percent level. The total one-halo 
correction P lh in real space is 7-10%. If we conservatively assume 
that the halo reconstruction algorithm incorrectly subtracts the real 



space one-halo term by 20%, then the systematic error at k = 
0.1 ft Mpc" 1 , ito.i, is allowed to be 2%. At k = 0.2 ft Mpc" 1 , 
the same error would translate to 5% in real space, though in red- 
shift space this term is mitigated. In Appendix IB4I we find that 
the shape difference of Phaio(k) and PLRo(k) is only 18% at 
k = 0.2 ft Mpc" 1 , and only 8% after accounting for the shot noise 
introduced by the fiber collision corrections. If we assume that our 
modeling and treatment of the one-halo contribution to the FOGs 
are accurate at the ~ 50% level, we can estimate a conservative 
error at k = 0.2 ft Mpc" 1 of 5%. Therefore, for all the modeling 
uncertainties besides central velocity dispersion that we have dis- 
cussed, Uo.i = 0.02 and «o.2 = 0.05 encompass the estimated un- 
certainties. These are our 'fiducial' nuisance function constraints. 

In Appendix IB2l we find that a large amount of central galaxy 
misidentification or central-halo velocity bias can reduce the am- 
plitude of Phaio(k,p) by a smoothly varying function of k at a 
level that exceeds these fiducial bounds on tto.i and Mo. 2- Our 
approach to mitigating the impact of uncertain central LRG pe- 
culiar velocities is twofold. First, for all of the analysis in Sec- 
tion [5] we adopt more conservative bounds for the nuisance func- 
tion: uo.i = 0.04 and 110.2 = 0.10, which nearly encompass the 
change in power spectrum shape in Fig. lBll for the extreme velocity 
dispersion model. Furthermore, we calibrate a second model from 
the mocks with extreme velocity dispersion, and in Appendix Iclwe 
determine the cosmological parameter constraints with this model 
to establish the level of remaining systematic uncertainty in our fi- 
nal results. 



5 COSMOLOGICAL CONSTRAINTS 

In this Section we explore the cosmological constraints derived 
from the power spectrum of the reconstructed halo density field, 
Phaio{k). We first consider constraints obtained from Phaio{k) 
alone, and then combine the L RG likelihood with W MAP5 and 
the Union Supernova Sample (Kowalski et al. 2008) to explore 
joint constraints in several cosmological models. Throughout, we 
make use of the COSMOMC package dLewis & Bridle! 120021) to 
compute cosmological constraints using the Markov Chain Monte- 
Carlo method. A stand-alone module to compute the Phaio(k) like- 
lihood is made publicly available^ 



5.1 Constraints from the halo power spectrum 

In this subsection we examine the cosmological constraints de- 
rived from the Phaio(k) alone and in combination with a prior on 
Q m h 2 from WMAP5. In the model Phaio(k,p), the scale factor 
a sc .i in Eqn.[8]is evaluated at z e ff = 0.313. For comparison with 
other works, we scale our constraint on Dy (0.313) using the fidu- 
cial distance-redshift relation, for which Dv{0-35)/Dy(z e ff) = 
1.106; the variation of this ratio with cosmo logical parameters 
is negligible. Following lEisenstein et al. I d2005t) . we consider two 
free parameters Q m h 2 and Dv{0.35). In this subsection we hold 
fl b h 2 = 0.02265, n s = 0.960, and <t 8 = 0.817 fixed at their val- 
ues in the fiducial cosmological model, and assume a flat ACDM 
model; in ^5.2| we relax these assumptions. 

For the 45 bandpowers satisfying 0.02 < k < 0.2 ft Mpc" 1 , 
X 2 is minimized when ZV(0.35) = 1396 and f2 m ft 2 = 0.136 
with best-fitting nuisance parameters ax = 0.160 and a-z = 



2 http://lambda.gsfc.nasa.gov/toolbox/lrgdr/ 
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Table 2. One-dimensional constraints from the LRG P; la ; (fc) likelihood, or in combination with the WMAP5 U m h 2 constraint U m h 2 = 0.1326 ± 0.0063 
('+ prior', below the line). We vary the kmax (units of h Mpc -1 ) included in the fit, the nuisance function constraints (fiducial vs weak F n uis)< velocity dis- 
persion in the model (fiducial vs. 'VD'), and whether the BAO features are included in the model (fiducial vs 'NW'). All constraints have assumed the ACDM 
relation between il m , Hq, and Dy, Q b h 2 = 0.02265, n s = 0.96, and <rg = 0.817. In the last column we show A0.35 = yj f! m -f/Q-EV(0.35) /0.35c 
lEisenstein et alj200^) . Models with weak F nu i a constraints or central galaxy velocity dispersion are discussed in Appendix|C] The kmax = 0.2 h Mpc -1 
constraints highlighted in bold are our main results, and the other cases are shown for comparison. 
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Figure 5. Constraints from the LRG DR7 Phalo{k) for a ACDM model 
with f2(,/i 2 = 0.02265 and n s = 0.960 fixed. The dotted contours show 
A\ 2 = 2.3 and 6.0 contours for the Phalo(k) fit to a no-wiggles model. 
The solid contours indicate A% 2 = 2.3, 6.0, and 9.3 contours for k max = 
0.2/i Mpc -1 and our fiducial Phalo(k, p) model. The three dashed lines 
show the best-fitting and ±ler values r s /Dv(0.35) = 0.1097 ± 0.0036 
from P09. 

-0.181: x 2 = 39.6 for 40 de srees of freedom. Thus the assumed 
model power spectrum and covariance matrix provide a reasonable 
fit to the observed spectrum. In a ACDM model, this point corre- 
sponds to h — 0.67 and fi m = 0.30. Fig [5] shows \ contours 
in the Q m h 2 — Dv'(0.35) parameter space, while Table [2] reports 
marginalized one-dimensional constraints for several combinations 
of these parameters. 

The information in Phaio(k) can be roughly divided into 
broad-shape information and information from the BAO scale. 



Since in this subsection n s is fixed, the shape information is the 
location of the turnover in the power spectrum set by matter- 
radiation equality, which constrains Q, m h 2 Dv', information from 
the BAO scale constrains r s /Dv- Here, r s is the sound hori- 
zo n at the baryon-drag ep och, which we evaluate using Eqn. 6 
of lEisenstein & Hul |l998). These two scales correspond to con- 
strai nts on feO°„ 93 and feO ~°' 37 respectively, in a ACDM cosmol- 
ogy dTegmark et alj|200of) . 

To isolate information from the power spectrum turnover 
and exclude that of the BAO scale, we alter our model so that 
fdam P (fc, p) = P n w(fc, p) in Eqn. [TOl The dashed lines in Fig. [5] 
show the constraints when using this 'no wiggles' model. Most 
of the available shape information comes from large scales with 
k < 0.1 feMpc -1 ; we demonstrate this in Table|2]by fitting the 
Pda,mp(k, p) = Pnw(k, p) model with the Sl m /i 2 prior to the data 
up to kmax = 0.1 /iMpc -1 and k m ax = 0.2/i Mpc -1 . The num- 
ber of independent modes is proportional to (k max — k min ); thus 
between k — 0.1 h Mpc -1 and k — 0.2 h Mpc -1 there are about 7 
times more modes than between k m i„ and 0. 1 h Mpc - . Neverthe- 
less, the constraint on Q m h 2 Dv{0-35) only improves by w 10% 
with the inclusion of modes between k m ax = 0.1 /iMpc -1 and 
kmax = 0.2 /iMpc -1 and does not shift appreciably. This also in- 
dicates that our modeling in the quasi-linear regime 0.1 < k < 
0.2 h Mpc -1 does not bias or substantially improve this constraint. 

If we reintroduce the BAO features in the model Phaio(k, p), 
then the degeneracy between (0.35) and Q m h 2 is partially bro- 
ken (solid contours in Fig. [5), and the constraints grow tighter 
as we include additional modes. This is understandable as the 
region 0.1 < k < 0.2 h Mpc -1 includes the location of 
the second BAO. The constraints on both r s /ZV(0.35) and 
Q m h 2 Dv (0.35) listed in Table [2] improve with k ma x- The mean 
value of Q m h 2 Dv (0.35) is consistent with what we find using the 
fdam P (fc, p) = -Pnw(fc, p) model with the WMAP5 Q, m h 2 prior, 
and does not shift substantially with increasing kmax - Because the 
BAO features break the degeneracy between Q m h 2 and _EV (0.35), 
the LRG Phaio(k) provides an independent constraint on Q m h 2 . 
For n s = 0.96, we find Qmh 2 = 0.141+°'qj 2 , which is consistent 
with the WMAP5 constraint, Q. m h 2 = 0.1326 ± 0.0063, but with 
a 70% larger error. 
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Fig. [5] shows that the LRG-only constraints derived with 
kmax = 0.2 ft Mpc -1 are consistent with the intersection of 
the power spectrum shape constraint (dotted lines) combined 
with constraints on r s /IV (0.35) from P09: the best-fitting and 
±lcr lines, 0.1097 ± 0.0036 are shown as dashed lines. Note 
that these are one parameter la errors. Table [2] shows excellent 
agreement for this quantity for the LRG-only constraints, with 



3 /£V(0.35) = 0.1097 



+0.0039 



for k ma x = 0.2 ft Mpc .This 



agreement reinforces the argument in Appendix IA2I that our ne- 
glect of the model dependence of the window function does not 
introduce significant bias in the £V(0.35) constraint. Moreover, 
this constraint does not change if we adopt very weak constraints 
on the nuisance function, \F nuis (0.1 ft Mpc -1 ) l/fa 2 , < 0.2 and 
l-Fn^fO^ftMpc- 1 )!^ < 0.5, or use the extreme central galaxy 
velocity dispersion model instead. We show in Appendix [C] that 
the largest known source of systematic uncertainty, the central 
galaxy velocity dispersion, impacts the cosmological parameter 
constraints at well below the statistical errors, and can be safely 
neglected for this analysis. We also demonstrate that our results are 
robust to the treatment of the nuisance parameters ci\ and 02- 

We estimate the significance of the detection of the BAO fea- 
ture as the difference between the best-fitting \ 2 for the fiducial and 
no wiggles models when ai, 0.2, and b 2 are chosen to minimize x 2 ; 
we find Axbao ~ 8-9- The resulting constraint on r s / Dv(0.35) 
is much tighter than is available from the shape information alone. 
To see this result, in Table[2]we combine the LRG Phaio(k) likeli- 
hood with a WMAP5 prior onf! m h 2 . The constraint from the shape 
alone, obtained by fitting the no wiggles model, gives a constraint 
on r s /Dv that is consistent with the constraint from the model in- 
cluding BAOs, but with a factor of ~ 2.3 larger errors. Finally, we 
note that P09 estimate the total BAO detection significance to be 
A\ 2 = 13.1; it is substantially larger than the value we find due to 
the inclusion of lower redshift galaxies from both the SDSS main 
sample and 2dFGRS. 

Finally, Table [2] also reports our constraint on A0.35 
jEisenstein et al. 2005): 



An 



H 2 



JV(0.35) 
0.35c 



(22) 



This parameter is tightly constrained by the Phaio(k) measurement 
and is independent of Ho- 



5.2 Dependence of LRG-only constraints on the cosmological 
model 

In Section lBTl the cosmological parameters Qbh 2 , n s , and as were 
fixed at their WMAP5 recommended values. For our purposes, r s 
changes negligibly as a function of flth since this parameter is 
so tightly constrained by CMB data. The parameters Q, m h 2 and n s 
both affect the linear power spectrum and are degenerate in shifting 
the contours along the constant r a /Dv direction, as illustrated in 
the upper panel of Figure [6] This degeneracy is well described as 
ttmh 2 (n a /0.96) 1 ' 2 = 0.141. 

In Figure [5] we have assumed the ACDM relation between 
fi m , ft, and Dy. This determines the scale at which to apply 
the non-linear corrections, which are at fixed k values in units of 
ft Mpc -1 . In the bottom panel of Figure |6]we show that this as- 
sumption is not restrictive. The dashed curve fixes ft = 0.7 and 
assumes no relation between ft and Dy. which also depends on Q,k 
and w. Varying as by ±0.1, which enters the HALOFIT calculation 
of the smooth component of the non-linear matter power spectrum 
in Eqn.[TT] changes the contours in Figure[5]negligibly. 
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Figure 6. Upper panel: Change in A\ 2 = 2.3 contour as n s is varied, 
with all other parameters as in Fig. [5] n s = 1.02 (dashed), n s = 0.96 
(solid), and n B = 0.90 (dotted). The degeneracy is well-described as 
Qmh 2 (n a /0.96) 1 ' 2 = 0.14. Lower panel: The impact of assuming a 
ACDM relation between f2 m , ft, and Dy (solid contours) compared with 
applying the non-linear corrections at ft = 0.7 and assuming no relation 
between Q m , h, and Dy (0.35) (dashed contours). As in Figf5]the lines 
show the constraints for constant r s /Dy (0.35) from P09. 



5.3 Combined constraints with WMAP5 and Union SN 

As probes of the redshift-distance relation, the three cosmologi- 
cal datasets we use in this Section are highly complementary for 
constraining the geometry of the universe and the equation of state 
of dark energy: WMAP5 effectively constrains the distance to the 
surface of last scatter and Q m h 2 , supernova data constrains angu- 
lar diameter distance ratios up to z ~ 1, and Phaio(k) sets joint 
constraints on r s /ZV(0.35) and f2 m ft 2 (n a /0.96). In Fig. [7] we 
show the intersection of these constraints for two models assum- 
ing a power law primordial power spectrum and no massive neu- 
trinos. The blue bands indicate the WMAP5 constraints and the 
gree n bands show constr aints using the Union Supernova Sam- 
ple dKowalski et al ] l2008h . For the P ha io(k), we show the con- 
straint on A0.35 (open bands), which has assumed n s = 0.96 and 
Qih 2 — 0.02265, and is independent of Hq. In the upper panel, we 



© 0000 RAS, MNRAS 000, 000-000 



Cosmology from the SDSS DR7 LRG Clustering 13 




Figure 7. WMAP5, Union supernova sample, and the LRG Phalo(k) 
Afj.35 constraint on the geometry of the universe. Upper panel: curva- 
ture varies and to = — 1 is fixed. The dashed line shows a flat universe, 
Q m + SI a = 1- Lower panel: w varies (assumed independent of redshift), 
and a flat universe is assumed. The dashed line indicates a cosmological 
constant, uu = — 1, WMAP5 and Union supernova contours are MCMC 
results, while for i\ a ; (fc), we approximate A\ 2 = 2.3 and A\ 2 = 6 
contours by showing A0.35 ± V2.3ovt 35 and Ay. 35 ± v / 6.0<t j 4 35 from 
the constraints in the top row ofTable[f| 

have assumed w = — 1 and allow curvature to vary. The three in- 
dependent constraints intersect near Q m = 0.3 and a flat universe 
(dashed line). In the lower panel, we assume flatness but allow w 
to vary; again the contours intersect near Q m — 0.3 and w = — 1, 
a cosmological constant. 

In this Section we combine these probes using the Markov 
Chain Monte Carlo (MCMC) method to obtain constraints on four 
cosmological models: a flat universe with a cosmological constant 
(ACDM), a ACDM universe with curvature (oACDM), a flat uni- 
verse with a dark energy component with constant equation of 
state w (wCDM), and a wCDM universe with curvature (owCDM). 
In each model we combine the constraints from Phaio{k) with 



the WMAP5 results jPunklev et~al]|2009l) . In the last model, we 
also present constraints in c ombination with both WMAP5 and 
the Union Supernova Sample dKowalski et al.ll2008h . Marginalized 
one-dimensional parameter constraints are presented in Table [3] 

The best-fitting ACDM fit to the WMAP5+LRG likelihoods 
is (Q m , Ob, Ov, Us, a 8 , h) = (0.291, 0.0474, 0.709, 0.960, 0.820, 
0.690) with best-fitting nuisance parameters ci\ — 0.172 and 
a 2 = -0.198. This model has Xlrg = 40.0 when fitting to 45 
bandpowers, and is shown with the data in Fig. [8] In this model 
adding the information from Phaio breaks the partial degeneracy 
between Q m and Ho in the WMAP5 data and reduces the uncer- 
tainties in each by a factor of ~ 1.6 compared to WMAP5 alone: 
Q m = 0.289 ± 0.019 and H = 69.4 ± 1.6 km s" 1 Mpc" 1 
(Q m = 0.258 ± 0.03 and H = 71.9±f% km s' 1 Mpc" 1 for 
WMAP5). The constraint on os also tightens by 30% because of 
the as — Q m h 2 partial degeneracy in the WMAP5 data. Note that 
since we marginalize over the galaxy bias, we have no constraint 
on as directly from the LRGs. 

In Fig. [9] we show the effect of opening the cosmological pa- 
rameter space to include curvature and a constant dark energy equa- 
tion of state w. Solid contours show the ACDM constraint in each 
panel for comparison. The dashes show WMAP5-only constraints. 
Without the ACDM assumption, WMAP5 cannot constrain Q m 
and Ho separately from Q m h 2 . In each of these models, the inclu- 
sion of the Phaio(k) information can break the degeneracy through 
the BAO constraint on r s /Dv- Table [3] shows that the cold dark 
matter density, Q c h 2 , constraint improves by ~ 15% compared 
to the WMAP5-only constraint (~ ±0.0063) due to the power 
spectrum shape information in the non-ACDM models. Moreover, 
the r s /ZV (0.35) constraint does not deviate substantially from the 
Phaio(k)+O.mh 2 prior constraint presented in Table[2] In the con- 
text of power-law initial conditions, Phaio(k) information does not 
improve constraints on the spectral index n s . 

Allowing curvature relaxes the constraints on Q m and Hq to 
the WMAP5-only ACDM errors on these parameters, while tightly 
constraining Q t ot = 1 - to 1.0114±°;°°^ (-0.027 < Q k < 
0.003 with 95% confidence). If instead we assume flatness but al- 
low the dark energy equation of state as an additional parameter w 
(assumed constant), w is constrained to —0.79 ± 0.15. Since the 
effective LRG sample redshift is z e ff = 0.313, allowing w to de- 
viate from —1 significantly degrades the 2 = constraints, Q m 
and Ho- 

When both fit and w vary, there remains a large degeneracy 
between Q m , Ho, and w. Curvature is still tightly constrained and 
consistent with flatness at the percent level: Qtot = 1.009 ± 0.012. 
Figure [10] demonstrates that supernovae can break the degeneracy 
in this model. The combination of all three data sets simultaneously 
constrains Qk within 0.009 and w to 11%, while still improving 
constraints on Q m and Ho compared with WMAP5 alone in the 
ACDM model. Allowing Qk 7^ and/or w 7^ —1 all act to increase 
Q m and decrease Ho compared with the ACDM model. The upper 
panel of FigureQJJshows that the ACDM model is only ~ la away 
from the best fit. The full set of constraints on all parameters is 
reported in Table [3] 

5.4 Additional constraints from the broad Phaio(k) shape 

For the models considered thus far, we have shown that gains in 
cosmological parameter constraints from adding constraints on the 
broad shape of Phaio(k) to WMAP5 results are moderate: ~ 15% 
improvement in Q c h 2 for all the models considered in Table [3] 
On the other hand, when the constraints on Qbh 2 and Q c h 2 from 
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Figure 8. Points with errors show our measurement of -Pf l(l ; (fc). We show \J Cu as error bars; recall that the points are positively correlated. We plot 
the best-fitting WMAP5+LRG ACDM model {Q, m ,Cl b , fi A , n s , erg, h) = (0.291, 0.0474, 0.709, 0.960, 0.820, 0.690) with best-fitting nuisance parameters 
a\ = 0.172 and = —0.198 (solid curve), for which \ 2 = 40.0; the dashed line shows the same model but with a\ = ai = 0, for which \ 2 = 43.3. 
The BAO inset shows the same data and model divided by a spline fit to the smooth component, P amoo th, as in Fig. 4 of P09. In Section |5~T1 we find the 
significance of the BAO detection in the Phalo(k) measurement is A\ 2 = 8.9. 



WMAP5 are used, our constraint on the BAO scale provides a much 
more precise determination of Dv at the effective redshift of the 
survey than the shape information alone. 

In more extended models than we have thus far considered, 
we may expect the additional shape information to allow tighter 
constraints. The cosmological parameters most closely constrained 
by the broad P(k) shape are those which affect the shape di- 
rectly or which affect parameters degenerate with the shape: these 
are expected to be the power spectrum spectral slope n s , its run- 
ning dn s /dlnk, neutrino mass m„, and the number of relativis- 
tic species N e ff. Thus far in our analysis, we have assumed 
dn s /d\n k — 0, m v = 0, and N e ; / = 3.04. 

One intuitively expects the measurement of Phaio(k) to im- 
prove constraints on the primordial power spectrum. In a ACDM 
model where both running of the spectral index and tensors are 
allowed, WMAP5 still places relatively tight constraints on the pri- 
mordial power spectrum: n s = 1.0871q'q73 ant ^ dlnn s /dlnk — 
—0.05 ± 0.03. The measurement reported in this paper probes 



at most A In k ~ 2 and covers a range corresponding to £ ~ 
300 — 3000; this range overlaps CMB measurements but extends to 
smaller scales. Over this fc-range and for this model, WMAP5 con- 
strains the P(k) shape to vary by ~ 8% from variations in the pri- 
mordial power spectrum. Due to the uncertainties in the relation be- 
tween the galaxy and underlying matter density fields, our nuisance 
parameters alone allow Phaio(k, p) to vary by up to 10— 14% over 
this region. Therefore we do not expect significant gains on n s or 
d In n s jd In k from our measurement. 

The effect of massive neutrinos in the CMB power spec- 
trum is to increase the height of the high I acoustic peaks: free 
streaming neutrinos smooth out perturbations, thus boosting acous- 
tic oscillations. In the matter power spectrum instead, neutrino 
free streaming gives a scale-dependent suppression of power on 
the scales that large scale str ucture measurements currently probe 
dLesgourgues & Pastotl 12006) . This makes these two observables 
highly complementary in constraining neutrino masses with cos- 
mology. 
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Figure 9. WMAP5+LRG constraints on Vl m h 2 , Q m , and Hq for ACDM (solid black contours), oACDM (shaded green contours), wCDM (shaded red con- 
tours), and owCDM (shaded blue contours) models. The first three panels show WMAP5-only constraints (dashed contours) and WMAP5+LRG constraints 
(colored contours) in the Q m h - Q m plane as the model is varied. In the lower right we show all constraints from WMAP5+LRG for all four models in the 
Om — h plane, which lie within the tight Q, m h' 2 ~ 0.133 WMAP5-only constraints. 



parameter 


ACDM 


oACDM 


wCDM 


owCDM 


owCDM+SN 




0.289 ± 0.019 


0.309 ± 0.025 


0.328 ± 0.037 


0.306 ± 0.050 


0.312 ±0.022 


#0 


69.4 ± 1.6 


66.0 ± 2.7 


64.3 ± 4.1 




65.6 ± 2.5 


£V(0.35) 


1349 ± 23 


1415 ± 49 


1398 ± 45 


1424 ± 49 


1418 ±49 


r s /Dy (0.35) 


0.1125 ±0.0023 


0.1084 ±0.0034 


0.1094 ±0.0032 


n in7H+ 0033 
0.1078_ 0034 


0.1081 ±0.0034 


Q k 




-0 0114+ 0076 




-0.009 ± 0.012 


-0.0109 ± 0.0088 


w 






-0.79 ±0.15 


-1.06 ±0.38 


-0.99 ±0.11 


n A 


0.711 ± 0.019 


0.703 ± 0.021 


0.672 ± 0.037 


70 n+0.057 

u. <uj_ Q5g 


0.699 ± 0.020 


Age (Gyr) 


13.73 ± 0.13 


14.25 ± 0.37 


13.87 ± 0.17 


14.27 ± 0.52 


14.24 ± 0.40 






1 0114+ ' 0077 

1 - ul -0 .0076 




1.009 ±0.012 


1.0109 ±0.0088 


100Of,/i 2 


2.272 ± 0.058 


2.274 ± 0.059 


9 9 Q o+0.062 
Z.Z30_q q 63 


9 9-Q+0.066 
z.zo_ Q65 


9 97f; +0.060 
z -^' D -0.059 


n c h 2 


H61+ 0039 
u.noi_ 0038 


0.1 1 10 ± 0.0052 


n , -,-,9+0.0056 
u - lllz -0.0057 


H03+ 0055 
u.riuo_ 0054 


ino+ 0051 
u.inu_ 0052 


T 


0.084 ± 0.016 


0.089 ± 0.017 


0.088 ± 0.017 


0.088 ± 0.017 


0.088 ± 0.017 


n s 


0.961 ± 0.013 


0.962 ± 0.014 


0.969 ± 0.015 


0.965 ±0.016 


0.964 ± 0.014 


ln(10 10 A)5) 


3.080lH 3 37 


3.068 ± 0.040 


o n71 +0.040 
o - u,1 -0.039 


3.064 ± 0.041 


3.068 ± 0.039 


C"8 


0.824 ± 0.025 


0.796 ± 0.032 


0.735 ± 0.073 


0.79 ±0.11 


n 7 go+ - 045 
u -' yu -0.046 



Table 3. Marginalized one-dimensional constraints (68%) for WMAP5+LRG for flat ACDM, ACDM with curvature (oACDM), flat wCDM (wCDM), wCDM 
with curvature (owCDM), and wCDM with curvature and including constraints from the Union Supernova sample. Here r is the optical depth to reionization, 
n B is the scalar spectral index, and A05 is the amplitude of curvature perturbations at k = 0.05/Mpc; these parameters are constrained directly by the CMB 
only. 
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Figure 10. For the owCDM model we compare the constraints from 
WMAP5+LRG (blue contours), WMAP5+SN (green contours), and 
WMAP5+LRG+S./V (red contours). In the upper panel the vertical line 
indicates a flat universe (Cj. = 0), and the horizontal line indicates 
a cosmological constant (w = —1). In the both panels we overplot 
the WMAP5+SN 68% contour (solid black) and WMAP5+LRG (dotted 
black) for ease of comparison. 



We start by comparing the constraints from 
WMAP5+P halo (k) and WMAP5+BAO (using the P09 BAO 
likelihood) in the ACDM model with three degenerate massive 
neutrino species. While WMAP5 alone finds ^2 m v < 1.3 eV 

with 95% confidence, WMAP5+Aai (fc) yields J2 m " < °- 62 
eV, which is a significant improvement over ^ m v < 0.73 eV 
(WMAP5+BAO). The upper panel of Fig. QT| compares the likeli- 
hood for m v for WMAP5 data alone (dashed) and in combination 
with P ha i (k). 

A change in the number of relativistic species in the early 
universe changes the epoch of matter-radiation equality and thus 
shifts the CMB acoustic peaks. The CMB constrains the redshift 
of matter-rad iation equality throu gh the ratio of the third to first 
peak heights dKomatsu et al.ll2009l) . If the effective number of rel- 



Figure 11. Phalo{k) improves constraints on neutrinos in the ACDM 
model through both the BAO scale and the broad power spectrum shape 
constraints. We show the one dimensional cumulative probability for 
WMAP alone (dashed) and WMAP+P haio (fc) (solid) for the neutrino mass 
(upper panel) and the one dimensional likelihood for the effective number 
of relativistic species N e f i (lower panel). 



ativistic species N e f / is allowed to va ry, this constraint de fines a 
degeneracy between Q, c h 2 and JV e // (Dunkley et al. 2009). Note 
that the physical quantity that is being constrained is the physi- 
cal energy density in relativistic particles. In the standard model 
this is given by photons and neutrinos but N v should really be 
considered an "effective" number of relativistic neutrino species: 
N e ff — 3.04 for standard neutrinos. Departures from this number 
can be interpreted also in terms of decay of dark matter particles, 
quintessence, exotic models, and additional hypothetical relativis- 
tic particles such as a light majoron or a sterile neutrino. 

In the ACDM model, which specifies a rigid relation between 
the angular diameter distance at last scattering measured by the 
CMB and low redshift distance scales, the degeneracy between 
N e f f and Q c h 2 can be broken by a low redshift distance constraint 
such as the BAO. However, N e f / will also impact the matter power 
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spectru m, which probes the hor izon size at matter-radiation equal- 
ity (e.g.. lEisenstein & Hul l 1998). Therefore, Phaio(k) is an excel- 
lent probe of N eff : WMAP5+P ha i (k) finds N eff = i.8±{%, 
while WMAP5+BAO yields N eff « 6.0 ± 2.5. For comparison, 
Koma tsu et al. I d2009h find N e ff = 4.4 ±1.5 when combining 
WMAP, BAO, supernovae, an d the Hubble Space Telescope key 
project (Freedman et al. 2001). The lower panel of Fig. [TT] com- 
pares the likelihood for N e f / for WMAP5 data alone with a prior 



N, 



cff 



10 (dashed) and in combination with Phaio{k)\ no prior 



on N e f f is needed in this case. 



6 COMPARISON WITH OTHER ANALYSES 

6.1 Comparison with previous galaxy clustering results 

There have been several previous an alyses of the clusterin g of 
the SDSS LRG spectroscopic sample. [Eisenstein et al.l j2005h use 
the correlation function of the DR3 SDSS LRG sample to de- 
rive constraints on Q. m h 2 = 0.133(n s /0.96)~ 1 ' 2 ± 0.011 and 
IV (0.35) = 1381 ± 64 Mpc, where we have adjusted their con- 
straints to match our assumed values of Qth 2 and n s ; recall that 
these constraints are not independent. Comparison with their Fig- 
ure 7 indicates that our model is slightly more than la away from 
their best fit. Our analysis prefers larger Q m h 2 and lower r s /Dv- 
In interpreting this comparison one should consider the differences 
in modeling and the fact that we have a factor of ~ 2 larger volume. 
Given this larger volume, naively we would expect an improvement 
on the constraints by a factor of ~ v2. Comparison with Table [2] 
shows that our LRG-only constraints on Q. m h 2 and Dy ha ve ap- 
proximately the same uncertainty as Eisenstein et al. (2005). This 
is partly because we conservatively increased our covariance ma- 
trix by a factor of 1.21 to account for the non-Gaussianity in the 
BAO contribution to the likelihood surface (see Section |2"l4l discus- 
sion). However, this increase will artificially weaken the constraint 
from the shape. Marginalization over the two nuisance parameters 
ai and ai to account for our uncertainty in the Phaio(k, p) as well 
as our conservative cut at k m i n also slightly weaken the constraint 
from the power spectrum s hape. 

iTegmark et al.l d2006h report cosmological constraints from a 
somewhat larger LRG sample (SDSS DR4) and combine their re- 
sults with WMAP3 data. To compare LRG-only c onstraints, we use 
th e value derived f r om th e ITegmark et alj d2006l) power spectrum 
in lSanchez&Colel d2008l) : Q m h = 0.173 ± 0.017 for n s = 1.0 
and h = 0.72. For a ACDM model scaled to n 3 = 1.0, our 



+0.012 



Restricting ou r 



LRG-only constraints yield Q, m h = 0.200_ 

analysis to k max = 0.1 /iMpc" 1 to match Teg mark et alj f2006). 
we find Q m h — 0.195 ± 0.013. Besides the increase in sam- 
ple volume, the discrepancy between these results could be due to 
differences in the FOG compression and the degeneracy between 
their nuisance parameter Q (see Eqn. Q} and cosmological param- 
et ers. A detai led comparison o f our modeling approaches is given 
in lReid et alj J2008I) . Note that lSanchez et all d2009h have also re- 
cently completed an analysis of the LRG correlation function, but 
they do not present a constraint from their shape measurement with 
which we can compare. 

Our results agree with analyses of photometric LRG samples. 
IPadmanabhan et alj d2007l) find n m = 0.30 ± 0.03 for h = 0.7 
and n, = 1. and lBlake etal] d2007» find Q m h = 0.195 ± 0.023 
for h = 0.75 and n 3 = 1. Our constraint is also consistent with de- 
terminations from other galaxy samples. For the 2dFGRS sample, 
ICole et al.1 {2005) find tt m h = 0.168 ± 0.016 for fixed n s = 1.0 



and h = 0.72; allowing a 10% Gaussian uncertainty in h yields 
Q m h = 0.174 ± 0.019, which is within la of our LRG-only con- 
straint. Our results are also i n good agreement with the SDSS main 
sample : ITegmark et aT] {2004b ) find Q m h = 0.201 ± 0.017, again 
with fixed n s = 1.0 and h = 0.72. 



6.2 Comparison with P09 

The P09 constraints overlap significantly with our analysis. 
We showed in Section 15.11 that our LRG-only constraint on 
r„/Dv (0.35) is in very good agreement with the determination 
in P09. When combined with the WMAP5 constraint on Q, m h 2 , 
our use of the shape information in Phaio(k) allows ~ 10% im- 
provement on r s /Dy (0.35). Moreover, the shape information pro- 
vides a tighter constraint on Q c h 2 . However, the P09 inclusion of 
SDSS main and 2dFGRS galaxies allows an additional constraint 
on r s /Dv (0.2), which generally makes the P09 constraints on fi m 
and Hq tighter. Our constraints on Q,k and w are comparable to P09. 
Across the models we have studied, WMAP5+Phaio(k) constraints 
yield lower values of Ho than the P09 results. This is driven by the 
P09 r s /Dv (0.2) constraint, which pulls the overall distance scale 
slightly lower compared to r s /Dy (0.35) alone, but does not signal 
any inconsistency between these analyses. 



6.3 Comparison with lRiess etaU d2009l) H 

iRiess et alj d2009ri recently released a new determination of the 
Hubble constant using a differential distance ladder: Ho — 74.2 ± 
3.6 km s _1 Mpc -1 . This value is consistent at the ~ la level 
with the WMAP5+P ha io(k) result for the ACDM model, H = 
69.4 ± 1.6 km s" 1 Mpc _1 . Tabled shows that if we allow fi fc / 
and/or w ^ — 1, the mean value of Ho decreases to ~ 64 — 67 
km s _1 Mpc -1 . Therefore, combining the lRiess et al.1 (2009) mea- 
surement with our constraints should reduce the uncertainties fur- 
ther and push the best-fit ting model closer t o ACDM. P09 present 
constraints including the IRiess et a[| {2009) Ho constraint for the 
owCDM model; the impact should be similar when using Phaio(k) 
rather than the P09 BAO constraints. 



7 CONCLUSIONS 

In this paper we have presented the power spectrum of the recon- 
structed halo density field derived from a sample of Luminous Red 
Galaxies (LRGs) from the Sloan Digital Sky Survey DR7. The size 
of LRG DR7 sample has sufficient statistical power that the de- 
tails of the relation between LRGs and the underlying linear den- 
sity field become important and need to be reliably modeled be- 
fore attempting a cosmolog ical interpretation of the data. Here, we 
have adopted the method of lReid et alj ( 2008), which applies a pre- 
processing step to the measured galaxy density field to reconstruct 
the halo density field before computing the halo power spectrum. 
On the scales of interest, this power spectrum has a more direct and 
robust connection to the underlying linear, real space power spec- 
trum than the power spectrum of the LRG galaxies themselves. 

We calibrate our method using Af-body simulations with vol- 
ume and resolution suitably tuned to trace the halo mass range rele- 
vant to LRGs, and provide several consistency checks between the 
observed and mock galaxy density fields to support our approach 
to model the LRG sample's clustering properties. In particular, we 
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demonstrate the validity of our modeling of the small-scale clus- 
tering and FOG features by matching the observed and mock cat- 
alogue higher-order statistics probed by the Counts-in-Cylinders 
group multiplicity function as well as the relative line of sight ve- 
locities between galaxies occupying the same halo. We discuss and 
quantify the sources of systematic error remaining in our model- 
ing. For the LRG sample, with nP ~ 1, both the shot noise sub- 
traction and the large velocity dispersions of their host haloes can 
introduce uncertainty. We identify the largest source of systematic 
uncertainty to be the velocity dispersion of central LRGs within 
their host haloes, and find its effects on cosmological parameters 
to be safely smaller than the size of the statistical errors. We are 
able to derive quantitative bounds on our model uncertainties and 
propagate these through the cosmological analysis by introducing 
nuisance parameters with tightly controlled allowed ranges, based 
on our understanding of the sources of non-linearity in the spec- 
trum. 

Based on our modeling of the LRG sample, we are able to 
extend our model for Phaio{k) to k = 0.2/iMpc - , increas- 
ing the number of available modes by a factor of ~ 8 over an 
analysis restricted to k ma x = 0.1 ft Mpc" 1 , as wa s the case in 
the SDSS team's DR4 analysis dTegmark et al]|2006l) . This allows 
us to simultaneously constrain the broadband shape of the un- 
derlying linear power spectrum and detect the BAO signal with 
A\ 2 = 8.9, though most of the shape information is confined to 
k < 0.1/iMpc" 1 . 

If we fix n s and Q,^h 2 , Phaio(k) alone constrains both 
Q m h 2 = 0.141±2loi2 and D v (0.35) = 1380±^. The agreement 
of our constraint on Qmh 2 at z e ff ~ 0.31 with the one derived 
from the CMB at z ~ 1000 provides a remarkable consistency 
check for the standard cosmological model. When Phaio{k) is 
combined with WMAP5, the error on Q c h 2 is reduced by ~ 15%, 
and the constraint on Dv"(0.35) allows us to place tight constraints 
on both Q m and Ho, as well as Qk or w. If we also include the 
Union Supernova Sample, all four parameters can be tightly con- 
strained: Q m = 0.312 ± 0.022, H = 65.6 ± 2.5 km s" 1 Mpc" 1 , 
fik = -0.0109 ± 0.008, and w = -0.99 ± 0.11, which is con- 
sistent with ACDM at the ~ la level. In fact, in the spirit of Oc- 
cam's razor, these constraints can be taken as evidence against both 
Ofe 7^ and w ^ — 1, since their values must conspire to match 
the observed angular diameter distance at recombination as well as 
£V(0.35); this can be seen from Fig. [7] 

Finally, we show that the shape information in Phaio(k) can 
improve constraints on both massive neutrinos and the number of 
relativistic species N e f f in a ACDM model. In combination with 
WMAP5 we find V m v < 0.62 eV at the 95% confidence level 
and N e ff = 4.8^]/*. These represent 16% (30%) improvements 
over using the WMAP5+BAO likelihood from P09. 

This paper represents a first attempt to analyse the LRG red- 
shift survey with a model that accounts for the non-linear galaxy 
bias and non-linear redshft space distortions introduced by the so- 
called one-halo term, and to propagate the uncertainty in the mod- 
eling through the cosmological constraints. We expect that the tech- 
nique introduced here to estimate the halo density field will be 
useful to furth er refinements such as reconstruction of the baryon 
acoustic peak dEisenstein et al. 2007a) and measurement of /3 from 
redshift space distortions. The modeling efforts presented in this 
paper are rather specific to the SDSS LRG sample. However, simi- 
lar techniques to probe the relation between the galaxy and under- 
lying matter density fields as well as to quantify its uncertainty will 
be required in the analysis of larger data sets from future galaxy 
surveys. 
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APPENDIX A: TESTING MODEL APPROXIMATIONS 

In this Appendix we present tests to demonstrate the validity of 
several assumptions of our model Phaio{k, p). 



Al Isotropy tests 

Both our Phaio(k,p) model (Eqn, |T6t and the a sc i approxima- 
tion (Eqn.[8} assume that the power spectrum modes are distributed 
isotropically with respect to the line of sight. We check this assump- 
tion in the SDSS DR7 LRG galaxy sample using pairs of galaxies 
separated by Ar m i„ = 15 h~ x Mpc to Ar m ax = 150 h~ 1 Mpc, 
binned into nine equal bins in Ar of width 15 ft -1 Mpc. We con- 
sider the two angles in the triangle defined by the observer and 
galaxy pair which give the angle between the galaxy pair separation 
vector and the local line of sight vector defined between the ob- 
server and one of the galaxies in the pair. These two angles will be 
equal in the limit of a pair with Ar <C max(xi, X2) where \i an d 
X2 are the distances to the two galaxies and Ar is the separation 
between them. We find (cos 2 (f>^ — 0.333 is —0.01 for the smallest 
separation bin (15 /i -1 Mpc < Ar < 30 /i -1 Mpc) and +0.005 
in the largest separation bin. Figure [AT] shows the full distribution 
versus | cos<j>\. The small increase for pairs perpendicular to the 
line of sight for the smallest separation bin is due to non-linear 
redshift space distortions (FOGs), inducing a potentially large sep- 
aration in redshift space between nearby pairs of galaxies in real 
space. The few percent deviations from isotropy will induce neg- 
ligible variations in the shape of the angle -averaged Ph aio(k,Tp), 
since the lower left panel of Fig ure 7 in lReid et al.U2008l) indicates 
only a ~ 5% change to the power spectrum shape between real and 
redshift space at k — 0.2 h Mpc . 



A2 Dy approximation 



As in Section [341 we use the approximation that pairs of galaxies 
contributing to Phaio(k) in the fc-range of interest are located at the 
same redshift to compute the effective survey redshift: 



Jzn 2 (z) 



m (z) dV 
b'-'(z) dz 



dz 



fn2(z) 



(Al) 



where n(z), b(z), and w(z) specify the average number den- 
sity, bias, and weight of the sample at redshift z as defined in 
IPercival et ail J2004) . We find z eff = 0.313, and use this redshift 
to evaluate a sc i in Eqn. [8] The effective redshift changes by only 
A 2 = 0.004 if one instead weights by the expected number of 
galaxies at redshift z. Given the distribution of pairs in the small 
separation limit (Eqn.ll8t we estimate the fractional bias remain- 
ing after the correction in Eqn.[8]is applied as 



8D\ 
Dy 



n 2 {^%dz 



J n 2 (z) 



w 2 (z) dV 
b 2 (z) dz 



dz 



-(A2) 



For a ACDM model, the fractional bias on the distance scale is < 
0.1% in the range Q. m = 0.2 — 0.4 and the rms change is < 1.2%. 
This additional variance about the peak is negligible for the BAO 




Figure Al. P(\ cos <f>\) vs cos (f> where </> is the angle between the galaxy 
pair separation vector and the line of sight defined by the observer and one 
of the galaxies in the pair (see text). The smallest separation bin (15 h 
Mpc < Ar < 30 hT 1 Mpc) shows the largest deviation from isotropy, with 
a ~ 5% preference for pairs perpendicular to the line of sight compared to 
along the line of sight due to FOGs. The larger separation bins extend to 
150 h~ 1 Mpc and are nearly isotropic, but with a few percent excess of 
pairs directly along the line of sight. 



scale ~ 100 ft -1 Mpc since this is much smaller than the damping 
scale o bao ~ 9 /i -1 Mpc. We find very similar results for the bias 
and rms damping if we instead integrate over the full distribution of 
isotropic pairs instead of using the Dy approximation in Eqn. lA2l 
Testing this approximation in more general models is more 
subtle, since Dv{z) depends on Ho, !J m , f2fc, and w, We in- 
stead do a consistency check: for f2 m /i 2 constrained by WMAP5, 
Dv{z e ff) constrained by WMAP5+/\ a ! (fc), 0. k = 0, and — 2 < 
w < —0.5, the maximum fractional bias is ~ 0.5%, and the max- 
imum rms change is 3.5%; a similar analysis for —0.025 < Qk < 
0.025 and w = — 1 shows much smaller deviations. We therefore 
conclude that in the range of models considered here, a single scale 
factor a sc l can accurately account for the effects of the model red- 
shift distance relation on the interpretation of the measured power 
spectrum. 



A3 Comparing P nw approximations 

In the models withou t massive neutrinos, we have used the 
lEisenstein & Hul ( 1998) formula (Eqn. 29) to compute P nw , which 
enters our model in Eqn. [10] However, for more general models 
such as those containing massive neutrinos or which vary the num- 
ber of relativistic species, it is more convenient to use a spline to 
obtain a smooth version of Pn n without BAO features. We fit a 
cubic b-spline to Pii n fc 1,5 in order to minimize the slope in the 
k region of interest. There are eight equally spaced nodes start- 
ing at k = 0.0175 Mpc -1 and ending at k = 0.262 Mpc -1 , and 
an additional node at k — 0.0007 Mpc -1 . Note we fix the loca- 
tion of the nodes in units of Mpc -1 since the linear power spec- 
trum is fixed in those units for fixed Q. m h 2 and Qth 2 . Fig. IA2I 
shows that the LRG-only likelihood surfaces computed with these 
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Figure A2. Comparison of the LRG-only likelihoo d surface comput e d with 
the analytic approximation of Eqn. 29 of lEisenstein & Hul jl998l) 

(solid, as in Fig. [5} compared with the result when using the b-spline fit 
described in Appendix I A3 1 (dashed) . 



two approximations agree well in the region preferred by WMAP5: 
Q m h 2 =0.133 ±0.0063. 



APPENDIX B: QUANTIFYING MODEL UNCERTAINTIES 
AND CHECKS FOR SYSTEMATICS: DETAILS 

In this Appendix we aim to quantify the sources of systematic un- 
certainty in the mode l Phaio(k). The model is calibrated on the 
mock catalogues of iReid et al] [2008). In Appendix IB 1 1 and Ap- 
pendix IB2I we present the detailed assumptions we have made to 
produce the mock catalogues from the TV-body simulation halo cat- 
alogues, and discuss the expected impact of these assumptions on 
the predicted relation between the reconstructed halo and matter 
density fields. Appendix |B3| through |B6| present consistency checks 
between the observed and mock catalogue LRG density fields that 
address the modeling uncertainties. In Section |4~3l the results of 
these tests are used to establish quantitative bounds on the nuisance 
parameters in Eqn. [T5] to be used in our cosmological parameter 
analysis. 



Bl Halo model parametrization 



In lReid et alj d2008l) we adopt the following para metrization for th e 
average number of LRGs in a halo of mass M dZheng et"ai]|2005h : 



(iVcen) = - 



(JV(AQ) = (iV c en> (1 + (AUt» 

Ifjgio M - lo Sl0 M ™ 



erf 



(N sat ) 



ClogM 
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(Bl) 
(B2) 
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For our adopted fiducial cosmological model, we find <ji os m ~ 
0.6—0.9 in order to match the amplitude of the observed large scale 
clustering of the LRGs; the exact parameter values used to generate 
the mock catalogues are given in lReid et al.l ( 2008). Since the scale 
dependence of halo bias varies with halo mass at the ~ 10% level at 



k — 0.15 h Mpc" 1 dSmith et all2007t) . changes in the distribution 
of LRGs with halo mass that preserve the large scale clustering 
amplitude could result in few percent changes in non-linear bias of 
the haloes traced by the LRGs. Changes in the distribution of halo 
biases traced by the LRGs could also alter the relation between the 
CiC and true group multiplicity function, which would introduce 
further uncertainty in the relation between the reconstructed and 
underlying halo density fields. 

B2 Distribution of mock galaxies within haloes 

In the mock catalogues of IReid et alj (2008) used to calibrate our 
model Phaio(k, p), we have assumed a sharp distinction between 
'central' and 'satellite' galaxies. The first or 'central' LRG in each 
halo is assumed to sit at the halo centre and move with the mean 
velocity of the halo dark ma tter; roughly 94% of the LRGs in ou r 
sample are central galaxies l lZheng et alj|2008l ; IReid et alj|2008l) . 
For the ~ 6% of LRGs that are 'satellites', we assume that they 
trace the phase space distribution of the halo dark matter, so that 
their positions and velocities are assigned to be that of a random 
dark matter particle in the halo. 

We do not evaluate the impact of errors in our assumed real 
space distribution of galaxies in their haloes on the fidelity of the 
halo density field reconstruction; the impact will be negligible in 
the case where there is a single LRG per halo. However, if the 
observed galaxies have a significantly different real space distri- 
bution in their haloes than we have assumed, the relationship be- 
tween the reconstructed halo density field and underlying matter 
density field will be different in the observed and mock galaxy 
catalogues. We test our assumed spatial distribution in Appendix 
|B3| by checking for consistency between the observed and mock 
catalogues for CiC group multiplicity functions, measured with 
two distinct sets of cylinder parameters. Furthermore, we can use 
Eqn.[20](where the measured CiC group multiplicity function spec- 
ifies {Ngai(N g ai — 1)}) as an upper limit on the error on the shot 
noise term due to differences between the model and observed re- 
constructed halo density fields. 

We consider two possible sources of deviation from our as- 
sumed galaxy distribution within haloes. The first is that on oc- 
casion an isolated LRG in our sample is not the 'central' galaxy 
in its halo, but a satellite galaxy, while the 'central' galaxy in that 
halo is not selected by our sample cuts. We call this situation 'cen- 
tral misidentification', and denote its probability fcen.err, assumed 
independent of halo mass for simplicity. The brightest LRGs are in- 
deed centrally concentrate^, wkh~ 80% of them within ~ 0.2r„i r 
of the X-ray peak dHo et al.|[2009T) . iLin & Mohj d2004l) similarly 
find that 80% of the BCGs in their X-ray selected cluster sample 
are within ~ 0.1r„i,-, and in the ~ 8% of cases where the BCG is 
outside 0.5r V i r , the second ranking galaxy in the group is within 
0.1r„i r . In some of these cases, both t he first and second brightes t 
galaxies would be identified as LRGs; Ivan den Bosch et alj (2007) 
showed that the luminosity difference between first and second 
brightest galaxies in massive groups is typically small. In this sit- 
uation there would be no error in our catalogues since we are not 
assigning luminosities to our mock LRGs. From these studies we 
would expect f cm ,err < 0.2 for the halo mass scales probed 



by these studies, M > 10 Mq, and it is reasonable to assume 
this holds at lower masses where there are fewer massive galaxies 
per halo. We therefore choose / ce „, e „ — 0.2 as our 'optimistic' 
value in the cases we consider in Fig. IB 1 1 U sing a galaxy group 



and c luster catalogue from SDSS dYang et al Jl 2007h . Tskibba et all 
find that the fraction of clusters in which the central galaxy 
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is fainter than the brightest satellite is « 30% in the mass range 



M ~ 10 1 



1O 14 M and w 40% for M 



10 14 - lO i& M . 
It is not clear what these results imply for the LRG galaxy sam- 
ple, but the parameter / ce „ ierr aims to encompass this case. We 
choose f C en,err = 0.4 as our 'conservative' estimate for the cases 
we consider in Fig. IB II 

The second situation we consider is the breakdown of our 
assumption that the central galaxy has no peculiar motion with 
respect to the mean velocity of the halo dark matter. Any offset 
with respect to the halo centre implie s that central galaxies are 
moving with resp ect to the halo centre ( van den B osc h etakll 20051 : 
Skibb a et allpred) . We call this situation central-halo velocity bias 
and parametrize the amplitude as b^ el = (Tc en /aoM> the ratio of 
the mean square veloci ty of the central galaxy to the halo dark mat- 
ter. Skibba et al. (prep) find b ve i ~ 0.1 once central misidentifica- 
tion has been accounted for. This small value is negligible for our 
purposes, so we set b ve i = in the 'optimi stic' and 'conser vative' 
cases we consider in Fig. IB II However, ICoziol et alj J2009h find 
bvei ~ 0.3 for brightest cluster members. This quantity is difficult 
to extract from observations, and it is not clear how the literature re- 
sults apply to the LRG sample because of the color-magnitude cuts 
defining the LRG selection. We set b ve i — 0.6 in the 'extreme' case 
we consider in Fig. IB II 

On the large scales of interest, the effect of nonzero f cen ,err 
or b ve i is to give the mock galaxies a velocity with respect to 
the halo centre. In Figure IB 1 1 we show the impact of nonzero 
central galaxy velocities on the recovered Phaio(k,p) for the 
three cases we described above. In the 'optimistic' case, we set 
{fce.n,err,b ve i) = (0.2,0); in the 'conservative' case, we set 
(fcen,err,b ve i) = (0.4,0); and in the 'extreme' case, we set 
{fcen.err, bvei) = (0.2, 0.6). To construct mock catalogues in each 
of these cases we leave the real space distribution of galaxies fixed. 
To mimic central misidentification, we replace the central galaxy's 
velocity with the velocity of a randomly selected dark matter par- 
ticle halo member. For central-halo velocity bias, we replace the 
central galaxy velocity with b ve iV ran , where v ran is the velocity of 
a randomly selected dark matter particle halo member. For compar- 
ison, we also outline both our fiducial nuisance function constraints 
(2% deviation at k = 0.1 /iMpc" 1 and 5% at k = 0.2 hMpc" 1 ) 
and the conservative nuisance function constraints (4% deviation at 
k = 0.1/iMpc" 1 and 10% at k = 0.2 h Mpc" 1 ) established in 
Section |4~3l The 'optimistic' case is well within the fiducial nui- 
sance constraints, and the 'conservative' case is well within the 
conservative nuisance constraints. The 'extreme' case, however, ex- 
ceeds the conservative nuisance constraints for k > 0.17 h Mpc - . 
In Section [5] we also evaluate the cosmological parameter con- 
straints when Phaio(k,p) is calibrated using the power spectrum 
of the 'extreme' velocity dispersion model in order to derive a limit 
on the systematic errors on our final results. 



B3 Comparison of mock and observed CiC group statistics 

In Table IB II we present CiC group multiplicity functions normal- 
ized by the number of galaxies per sample for two sets of cylinder 
parameters: Ar± < 0.8 h^ 1 Mpc, Av p = 1800 km s _1 (these 
are the parameters used to define our CiC groups and reconstucted 
halo density field for Phaio(k, Pad)) and Ar± ^ 1.2 h' 1 Mpc, 



k h Mpc -1 

Figure Bl. We compare three models including central galaxy velocity dis- 
persion to our fiducial model with no central galaxy velocity dispersion 
(fcen,err = b ve i = 0) by showing the ratio of Phalo(k, Pad) for the 
models. The dashed curve has fcen.err = 0.2, b ve i = ('optimistic'); 
the long dashed curve has f ce n,err = 0.4, b ve [ = ('conservative'); the 
dash-dot curve has / ce ?i,err = 0.2, b ve i = 0.6 ('extreme'). The straight 
lines show our fiducial (solid) and conservative (dotted) nuisance parameter 
constraints determined in Section|4~3l 



The observed groups contain 2158 LRGs that were as- 
signed redshifts by th e fiber collision correction. According to 
iReid & SpergeJ (2009), where colors are used as a redshift indica- 
tor, up to ~ 36% of these may be erroneous assignments; correct- 
ing this would remove ~ 780 galaxies from the observed groups. 
We find 6.2% of the observed galaxies are 'satellite' galaxies us- 
ing the reconstructed haloes, or 5.5% if we apply a correction for 
erroneous fiber collision assignments, while our mock catalogues 
have 5.9%. The structures of the multiplicity functions are gener- 
ally similar. Since our mock catalogues we re designed to mat ch this 
measurement but for LRGs selected as in Zeha vi et alj d2005h , the 
level of agreement is as expected. We verify that the agreement ex- 
tends to the multiplicity function when we adjust the group finding 
parameter Ar± to be 50% larger. Accounting for the possible con- 
tamination from fiber collision corrections, which is likely to man- 
ifest mostly at that in general the observed dis- 
tribution is smaller than in the mock catalogues at all multiplicities 
and for both values of Ar± . This result may be understood as one 
or more of three possibilities: the mocks having too many satellites 
altogether, different amounts of contamination from interlopers due 
to errors in the small-scale two-halo redshift space correlation func- 
tion, or a tighter distribution of satellite galaxies about the central 
one in the mocks. An error of the first kind would result in no error 
in the reconstructed density field; errors of the other kinds would re- 
sult in small changes to the effective shot noise or FOG features in 
the density field. The last line in Table |BT] shows that the difference 
in the effective one halo term derived from the mock and observed 
catalogues is < 2% of the total shot noise correction. Since the dif- 



Av p = 1800 km s 1 . The second CiC multiplicity function com- ference between P lh n ga i measured at Ar± ^ 0.8 h 1 Mpc and 



puted with larger Ar± is used to demonstrate consistency between 
the mock and observed catalogues. If the observed satellite galaxies 
were significantly less concentrated than in our mock catalogues, 
we would detect these galaxies when Ar± is increased. 



Ar± ^ 1.2 h" 1 Mpc is less for the observed catalogues compared 
with the mocks, we cannot be missing significant contributions to 
P lh due to a less concentrated distribution of the satellite galaxies 
in the observed haloes compared with the simulated ones; rather, 
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Table Bl. The observed and mock catalogue CiC group multiplicity func- 
tions of groups with n group galaxies for our fiducial group finding param- 
eters Arj_ max = 0.8 h" 1 Mpc, Av p = 1800 km s —1 and for a bigger 
^ r l,moi = 1-2 h^ 1 Mpc. The final row shows the ratio of the one-halo 
shot noise f^flG fEqn. 1201 to the standard shot noise correction l/n^jia- 

the increase in the number of groups comes from the increase in 
contamination from galaxies residing in nearby haloes. Our final 
conservative nuisance parameter bounds, discussed in Section l4~3l 
allow an error of the order of 40% in the one-halo shot noise sub- 
traction. Also note that because the maximum line of sight sep- 
aration (Av p = 1800 km s" 1 or ~ 20 h~ Mpc) is so large, 
the model CiC multiplicity functions are nearly identical when we 
consider the model with 'extreme' central galaxy velocity disper- 
sion. Finally, adding some spatial dispersion of the central galaxies 
would slightly reduce the number of CiC groups for an otherwise 
fixed catalogue; this may bring the models and observations into 
even closer agreement. 

We compute the line of sight separation of galaxies in the same 
CiC group as a probe of the accuracy of our model galaxy veloci- 
ties at the high halo mass end, where there is more than one LRG 
per halo. The comparison is complicated by the presence of fiber 
collision corrected galaxies, since their redshifts are artificially set 
to that of another galaxy in their group. We discard all such groups, 
and discard an equal fraction at each n group in our mock sam- 
ple. The resulting distributions are shown in Fig. IB2I The fiducial 
mocks with no central galaxy velocity dispersion appear to fit the 
data better, though neither matches the observed sharpness of the 
rise at small separations. Note that the fiducial mock catalogues 
with no velocity dispersion are determined only by the observed 
Ncic (n g roup)', no free parameters have been adjusted to match 
the observed velocity distribution. This comparison indicates that 
the residual FOG features in the reconstructed observed and mock 
halo density fields will be in satisfactory agreement. 

B4 Comparison of Phaio(k) and Plrg (k) 

In this subsection we examine the difference between the observed 
redshift space monopole spectrum PLRo(k) (no density field pre- 
processing of FOG features) and the power spectrum of the recon- 
structed halo density field, Phaio(k), and compare with our mock 
galaxy catalogues. This comparison provides an additional consis- 
tency check between the mock and observed LRG catalogues, and 
quantifies the effect of the halo density field reconstruction step on 
the measured power spectrum shape. 
We consider 

AP(k) = P L R G (k) - b 2 re iP hato (k) (B4) 

where b re i is a constant that parametrizes the enhancement of the 
overall bias when satellite galaxies are included, since they oc- 
cupy the most highly biased regions. In real space on large scales, 
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Figure B2. Solid line with error bars is the observed probability that a 
galaxy has a member of its CiC group with a separation Ar^os along 
the line of sight for pairs of galaxies identified as pairs by the CiC crite- 
ria, once all groups containing a fiber collision galaxy are removed. Error 
bars indicate fractional errors of 1 / '%/iV (A(r lo s) . giving a sense of the 
Poisson level of uncertainty in the measurement without considering the 
contribution from cosmic variance. The dashed line is the expected distri- 
bution for our model with no central galaxy velocity dispersion, and the 
dot-dashed line is for the model with central galaxy velocity dispersion. 
Note that Ar^os = 1 h Mpc corresponds to An 115 km s — 1 for 
the redshift distribution of our sample. 

AP(k) would be a simple shot noise, but in redshift space we ex- 
pect the detailed AP(fc) to result from the transfer of power be- 
tween scales caused by the FOGs, making AP(k) dependent on 
the underlying power spectrum shape. We will ignore this possi- 
ble < 10% level modification to the expected AP(k) since we 
have demonstrated good agreement between the shape of the mock 
and observed halo power spectra. The lower short dashed curve in 
Fig. IB3I shows the predicted AP(fc) from our mock catalogues, 
the upper short dashed curve shows the predicted AP(fc) scaled 
by a factor of 1.5, and the solid curve shows AP(k) for the ob- 
served spectra. The observed AP(k) is clearly shallower than the 
predicted shape. 

A crucial difference between the observed and mock LRG 
density fields is the application of fiber collision corrections dis- 
cussed in Section |2~T1 in the observed density field. 2158 galaxies 
without spectra were added to the LRG sample and assigned the 
redshift of the nearest LRG, while the CiC group multiplicity re- 
sults in Table IbTI indicate that 6857 galaxies are 'satellite' galax- 
ies. Fir st, since ~ 36% of th e fiber-collision corrections are erro- 
neous dReid&Spergel 2009), we expect an additional shot noise 
of ~ 125(ft _1 Mpc) 3 from these galaxies, which are not repre- 
sented in our mock catalogues. Second, the fiber collision corrected 
galaxies that are physically associated with a neighboring LRG will 
change the distribution of AP(k) relative to the mocks because 
their line of sight separation from the neighboring galaxy has been 
eliminated. The long dashed curves in Fig. IB3I shows that we can 
match the observed AP(k) as a sum of the mock catalogue AP(fc) 
and a shot noise of 200(/i~ 1 Mpc) 3 . The AP(fc) for the observed 
spectra is consistent with a constant power for k < 0.2 ft Mpc -1 
and amounts to a significant difference between the two spectra: 



,obs(™) N bigjmock (n) 
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Figure B3. The solid curve is the difference between the observed spectra 
Phalo(k) and Plrg (k), the lower short-dashed curve is the predicted dif- 
ference from our simulated catalogues, and the upper short-dashed curve 
is the same curve but scaled by a factor of 1.5. The scale dependence of 
AP(fc) is smaller for the observed spectra than for the simulation results. 
Furthermore, there is some uncertainty in the appropriate value of b re i, 
which changes the shape of AP(fc). However, at high k, the prediction is 
robust to changes in b re ; since P(k) is small. The long-dashed curves show 
AP moc j.(fc) + 200(/i — 1 Mpc) 3 for several values of b re i. This demon- 
strates that the difference between Phalo(k) and Plrg(^) i s consistent 
with the difference measured in the simulated catalogues if the excess shot 
noise from fiber collisions is accounted for. Moreover, the difference be- 
tween the observed halo and LRG spectra is large compared with the statis- 
tical errors on the bandpowers. 



~ 8% at k = 0.1 ft Mpc" 1 and ~ 18% at k = 0.2/iMpc^ 1 . 
Therefore, differences in the preprocessing of the LRG density field 
can lead to changes in P(k) much larger than the statistical errors 
on the measurements, which could then be propagated to errors in 
the derived cosmological parameters. Note that the reconstructed 
halo density field is basically unaffected by errors in the close-pair 
fiber collision correction applied to the data, since these galaxies 
are all assigned to haloes containing other LRGs already. 

In summary, the difference between Phaio(k) and PLRc(k) 
can be understood once we account for the effects of fiber colli- 
sions, and the model predictions Phaio{k, p) are robust to any un- 
certainty associated with these effects. 



B5 The effect of luminosity- weighting on Phaio (k) 

A further subtle difference between the mock and observed halo 
power spectrum is that the mock catalogues were evaluated using 
a redshift snapshot with constant ulrg, and luminosities were not 
assigned to the mock LRGs; each reconstructed halo is weighted 
equally when computing the overdensity field. To verify that the lu- 
minosity weighting used to compute the Phaloik) does not signif- 
icantly alter the relative amplitude of the shot noise to total power 
compared with our mock catalogues, we recompute Phaio(k) from 
the da ta with luminosity-independent weights from lFeldman et all 
dl994h : 



Figure B4. Ratio of the power spectra c omputed using the weights in 
Eqn. |B6| to the standard Percival et al. 12004) weighting scheme after reseat- 
ing the overall normalization. We also overplot Pu„(k) / P n w(k) for our 
fiducial model to demonstrate no correlation between the small shifts in 
the measured power spectrum and expected BAO feature. Errors show the 
fractional errors on Phalo(k), \fCTi / Phalo{k) . 



w(r, L) 



(B6) 



b(L) = 1 



(B5) 



1 + PoUlrg 

where P a = 10000(/i _1 Mpc) 3 . Fig. |B4] shows the ratio of 
the observed spectra with our fiducial weights compared with the 
luminosity-independent weights. The good agreement even at large 
k where the power is small indicates there is no significant differ- 
ence from the shot noise subtraction between these two weightings; 
we find no statistically significant change in the power spectrum 
shape. Moreover, the change in the windowed theory power spec- 
trum due to the change in weights is negligible (< 0.1%), indi- 
cating that the window function will not be sensitive to the partic- 
ular weighting choices of Section [2] for reconstructed haloes con- 
taining more than one galaxy. Whil e the luminosity-weig hting is 
critical for the SDSS main sample dTegmark et ai]|2004bh . Fig. □ 
shows that the LRGs are close to volume-limited over much of the 
redshift range of the sa mple; it is therefore u nsurprising that the 
lFeldmanetal] ( ll994h and lPercival et al.l ( 120040 weighting schemes 
produce nearly identical power spectra for the LRG sample. 



B6 Checking the cosmological dependence of the model 

Our model uses HALOFIT to describe the cosmological parameter 
dependence of the non-linearity in the matter power spectrum, and 
is calibrated from TV-body simulations at the fiducial cosmology 
(Eqn. 1 13b- Below k = 0.1 /iMpc -1 , the dark matter power spec- 
trum is linear at the 1% level, apart from the BAO damping, and 
it is only ~ 15% larger than the linear one at k = 0.2 ftMpc" . 
Using the publicly available WMAP5 ACDM MCMC chain, we 
find Phaioflt(fc) / Pun(k) changes by ~ ±2% for k ^ 0.2 in the 
space of cosmologies allowed by the WMAP5 data alone; the error 
on this small correction will therefore be well below 1%. Therefore 
we expect the model of the non-linear matter power spectrum to 
be accurate at the < 1% level at k = 0.1 h Mpc -1 and ~ 1% at 
k = 0.2/iMpc" 1 . 
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Figure B5. Agreement between Phalo.w MAP?,{k) I P n w{k, pwMAP'i) 
measured from the catalogues in Reid & Spergel 1 2009) based on an N- 
body simulation z = 0.2 snapshot with WMAP3 cosmological parameters 
(points with error bars) vs. the model prediction from Eqn ll6l at Znbar = 
0.235. 

We use the LRG catalogues from iReid & Spergeil 120091) 
evaluated at the WMAP3 preferred cosmological parameters 
(fl m ,n b ,n A ,n s ,(T8,h) = (0.26,0.044,0.74,0.95,0.77,0.72) 
at z — 0.2 with Ltox = 1 h~ l Gpc to test the cosmological de- 
pendence of our model Phaio(k,p) in Eqn. 1131 We plot a mock 
catalogue power spectrum Phaio,WMAPsik) / P nw {k,p W MAP3.) 
against our Phaio(k,p) model predictions for a NEAR subsam- 
ple in Fig. |B5| to demonstrate the agreement both in the BAO fea- 
tures and overall shape of the deviation out to k = 0.55 h Mpc - . 
X 2 = 96.6 for 86 DOF (k 0.55) and x 2 = 29.1 for 31 DOF 
(k 0.2). This provides further evidence that the cosmological 
dependence of our model Phaio(k,p) is sufficiently accurate for 
the SDSS DR7 data, which probe a somewhat smaller volume. 



APPENDIX C: EFFECTS OF CENTRAL GALAXY 
VELOCITY DISPERSION AND NUISANCE 
PARAMETERS 

In Section|4]we established that the largest remaining known source 
of systematic uncertainty is the central galaxy velocity dispersion. 
To test the impact of this uncertainty on the cosmological con- 
straints, we reevaluate the Phaio{k) likelihood surface using the 
'extreme' velocity dispersion model in Appendix IB 2 1 to calibrate 
the model Phaio{k,p). The maximum likelihood points for the 
fiducial, no velocity dispersion model (cross) and the 'extreme' 
velocity dispersion model ('X') are shown in the upper left panel 
of Fig. |C1| The systematic shift in the contours between the zero 
and extreme central velocity dispersion model is small compared 
to the width of the Ax 2 = 2.3 constraint (dotted curve). When 
we marginalize over nuisance parameters 6 2 , a%, and a?, A\ 2 be- 
tween the maximum likelihood model values for the zero and ex- 
treme velocity dispersion models is ~ 0.3. If one instead adopts 
the Oi, a2, and 6 2 values which minimize x 2 , the shift decreases 
to Ax 2 ~ 0.1; the difference is because the preferred nuisance pa- 
rameters a± and a,2 in the no velocity dispersion model are closer 



to the boundary of the allowed values. These Ax 2 values are ap- 
proximately the same when considering a fit to the model with 
or without BAO wiggles. This shift is small compared to the sta- 
tistical errors, and since the velocity dispersion model considered 
is extreme compared with the availab le estimates in the literature 
( SkibbaetalJlp re pT: ICoziol et al.| [2009). we can safely neglect this 
systematic uncertainty in the present analysis. 

Within our fiducial nuisance parameter bounds and us- 
ing our fiducial model with no central galaxy velocity dis- 
persion, we have verified that the effect of the nuisance pa- 
rameters in Eqn. Q3] is small on the Phaio{k) cosmologi- 
cal parameter constraints. The preferred nuisance parameters 
are off-center in the allowed ai — a-x space, although not 
at the boundary: (F nuis (0.1 /i Mpc _1 )/&o) = 0.016 and 

((F„ tti8 (0.1/iMpc- 1 )-ii , n«i a (0.2ftMpc- :L ))/^) = 0.060, 
where we have computed a likelihood-weighted average over the 
DR7-only constraints. The upper right panel of Fig. IC1I shows 
Ax 2 = 2.3, 6.0, and 9.3 contours where ai — a.2 = and b 2 
is varied to minimize x 2 (dashed contours) compared to our fidu- 
cial marginalization over b 2 , a\ and ai (solid contours). Allowing 
nuisance parameters to account for our imperfect modeling induces 
both a small shift and widening of the likelihood surface. The dif- 
ference in the contours is negligible when x 2 is evaluated instead at 
the values a\ and a 2 that minimize x 2 ■ Therefore the hard boundary 
we impose in a\ — a 2 space does not seriously affect the likelihood 
contours, and ai and a,2 are not strongly degenerate with the cos- 
mological parameters constrained by Phaio{k) when ai and a 2 are 
tightly constrained by the arguments in Section l4~3l 

However, when one substantially relaxes the constraints on 
the nuisance function, the constraints from the power spectrum 
shape degrade. The lower right panel of Fig. IC1I shows how the 
X 2 = 2.3, 6.0, and 9.3 constraints relax when a\ and 02 are cho- 
sen to minimize x 2 such that F nU i B (k — 0.1 h Mpc _1 )/6o < 0.2 
and F nuis (k = 0.2 ftMpc"')/6o < 0.5. While the constraints on 
r s / Dv(0.35) are unchanged, the shape information is degraded. 
The effects of sc ale dependent halo bias are well below these al- 
lowed deviations dSmith et alj|2007l) . and we have argued that our 
reconstruction of the halo density field should leave much smaller 
uncertainties as well. The dashed contours in the lower right panel 
of Fig. IC1I show a further broadening of the constraints when ai 
and a,2 are varied without restriction to minimize x 2 - For compar- 
ison with the adopted nuisance restrictions, the bottom right panel 
of Fig. IC II also shows the regions where the best-fitting nuisance 
parameters satisfy \F nuis (k = 0.1 h Mpc" 1 )!/^ < 0.04 (solid 
lines) and \F nuls (k = 0.2 h Mpc _1 )|/&o < 0.1 (dashed lines). 
The width of this region is smaller than the statistical errors de- 
rived from the shape constraint, which are shown in the upper left 
panel. Consequently, it is unsurprising that our marginalized like- 
lihood contours with the fiducial nuisance restrictions deviate only 
slightly from the contours where a\ — a2 = 0. Finally we note that 
for the models with and without velocity dispersion, the likelihood- 
weighted best-fitting nuisance functions have small deviations from 
one at k — 0.1 h Mpc -1 (< 2%), the region containing most of 
the shape information. The two models differ in the quasi-linear 
regime: ( (F nula (0.1 h Mpc" 1 ) - F nuia (0.2 h Mpc" 1 )) /6§) = 
—0.033 for the velocity dispersion model and 0.060 without veloc- 
ity dispersion. However, we cannot distinguish between velocity 
dispersion and other modeling uncertainties to explain the shape of 
the nuisance function preferred by the data. Moreover, using the 
velocity dispersion model does not improve the overall x 2 of the 
fit. 
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Figure CI. Effects of velocity dispersion and nuisance parameters on constraints from the LRG DR7 Phalo(k) for a ACDM model. In each panel we hold 
f!(,/i 2 = 0.02265 and n s = 0.960 fixed. Upper left panel: A\ 2 = 2.3 and 6.0 contours for the Phalo{k) fit to a no-wiggles model with no central velocity 
dispersion (solid) and extreme velocity dispersion (dashed). The A\ 2 = 2.3 for the fiducial model with BAO features is shown for comparison by the dotted 
line. The cross shows the maximum likelihood point for our fiducial model, while the 'X' shows it for the extreme velocity dispersion model. The solid line 
indicates r s / Dy (0.35) = 0.1097, demonstrating that adopting the velocity dispersion model shifts the likelihood surface along constant r s / Dy (0.35). 
Upper right panel: A\ 2 = 2.3,6.0, and 9.3 contours. The solid contours use our fiducial marginalization over b 2 , oi, and <J2 (as in Fig. |5J, while in the 
dotted contours fix a± = 02 = and b 2 to the value which minimizes \ 2 ■ Lower left panel: The solid contours as in Fig. [5] while the dashed contours take the 
minimum x 2 value for which \F(k = 0.1 fcMpc _1 )|/f;Q < 0.2 and \F(k = 0.2 h Mpc _1 )|/feQ < 0.5. Lower right panel: The solid contours as in Fig.[5] 
while the dashed contours minimize \ 2 with no restrictions on a\ and 02. For comparison with the fiducial nuisance restrictions, the solid lines enclose the 
region where for the best-fitting x 2 , \F(k = 0.1 h Mpc~ 1 )|/6§ < 0.04 and the dashed lines enclose \F(k = 0.2 hMpc _1 )|/fe^ < 0.1. 



We conclude that, for this data set, the statistical errors are 
comfortably larger than the errors from modeling uncertainties. 
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